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ABSTRACT 

The Spitzer Extended Deep Survey (SEDS) is a very deep infrared survey within five well-known extragalactic 
science fields: the UKIDSS Ultra-Deep Survey, the Extended Chandra Deep Field South, COSMOS, the Hubble 
Deep Field North, and the Extended Groth Strip. SEDS covers a total area of 1.46 deg 2 to a depth of 26 AB mag 
(3a) in both of the warm Infrared Array Camera (IRAC) bands at 3.6 and 4.5 fim. Because of its uniform depth of 
coverage in so many widely-separated fields, SEDS is subject to roughly 25% smaller errors due to cosmic variance 
than a single-field survey of the same size. SEDS was designed to detect and characterize galaxies from intermediate 
to high redshifts ( z = 2-7) with a built-in means of assessing the impact of cosmic variance on the individual fields. 

Because the full SEDS depth was accumulated in at least three separate visits to each field, typically with six-month 
intervals between visits, SEDS also furnishes an opportunity to assess the infrared variability of faint objects. This 
paper describes the SEDS survey design, processing, and publicly-available data products. Deep IRAC counts for 
the more than 300,000 galaxies detected by SEDS are consistent with models based on known galaxy populations. 

Discrete IRAC sources contribute 5.6 zb 1.0 and 4.4 zb 0.8 nW m -2 sr -1 at 3.6 and 4.5 ^m to the diffuse cosmic 
infrared background (CIB). IRAC sources cannot contribute more than half of the total CIB flux estimated from 
DIRBE data. Barring an unexpected error in the DIRBE flux estimates, half the CIB flux must therefore come from 
a diffuse component. 

Key words: galaxies: high-redshift - infrared: galaxies - surveys 
Online-only material: color figures, machine-readable tables 
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1. INTRODUCTION 

The Spitzer Extended Deep Survey (SEDS) is an unbiased 
deep survey carried out by the Spitzer Space Telescope (Werner 
et al. 2004). The survey consists of integrations to 12 hr 
per pointing at 3.6 and 4.5 /x m with Spitzer ’ s Infrared Ar- 
ray Camera (IRAC; Fazio et al. 2004b) in five well-studied 
fields (E-GOODS-S. E-GOODS-N, EGS, UDS, and COSMOS / 
UltraVISTA) with a goal of covering a total of one square de- 
gree. To study variability, each SEDS field was visited three 
times, obtaining four hours/pointing on each visit. Visits were 
typically separated by six months. SEDS required a total of 
2103 hr, making it the single largest Exploration Science Pro- 
gram carried out by Spitzer to date. 

SEDS was designed from the outset to carry out a census 
of stellar mass and black holes as a function of cosmic time 
reaching back to the era of reionization. SEDS was also designed 
to measure galaxy clustering over a wide redshift range to 
provide a critical link between galaxies and their dark matter 
halos. These observations will be critical for testing galaxy 
formation models in the early universe. These outcomes are 
made possible by IRAC’s unique ability to sample the rest-frame 
visible light of galaxies and active galactic nuclei (AGNs) at high 
redshift. The resulting broad-band spectral energy distributions 
can be combined with stellar population models to give stellar 
masses and ages. Specifically, the SEDS design goals were to 
detect galaxies down to ~5 x 10 9 M G at z = 6. This sensitivity 
is necessary to robustly measure M* at that redshift, i.e., to 
detect the galaxies that dominate the global stellar mass density. 
Based on the cosmological hydrodynamic simulations of Dave 
et al. (2006), specifically their favored momentum-driven winds 
model, this is achieved at 26 AB mag in the 3.6 pm band. The 
SEDS depth was set accordingly. The SEDS design also enables 
several additional studies, including (1) galaxy evolution in the 
redshift range z = 1-6, (2) obscured AGN growth at high 
redshift, (3) AGN variability, and (4) measurement of the cosmic 
infrared background (CIB) spatial fluctuations. SEDS explores 
the universe at a depth never before achieved over such a wide 
area at mid-infrared wavelengths. The depth and area coverage 
achieved by SEDS is compared to that from other major Spitzer/ 
IRAC surveys in Figure 1. 

This paper is organized as follows: Section 2 describes the 
five SEDS fields with particular attention to prior IRAC imaging 
incorporated into the final SEDS data products. Section 3 
discusses the details of the SEDS observing strategy and data 
reduction, and Section 4 describes the source identification, 
photometry, and validation. Section 5 describes the SEDS 
catalogs. In Section 6 we present preliminary results, including 
deep mid-IR number counts and the infrared color distribution 
of IRAC-detected galaxies. Section 7 summarizes these results. 
All magnitudes are stated in the AB system. 

2. THE FIVE SEDS FIELDS 

Because a key driver of SEDS science is the need to derive 
stellar mass estimates and galaxy stellar mass functions, it was 
vital to select fields with substantial quantities of deep broad- 
band photometry. Of special importance are the availability of 
near-IR and visible imaging deep enough to match the IRAC 
observations reported here. An exhaustive search of potential 
target fields led to a choice of five widely separated premiere 
deep-sky survey regions. These are the 30' x 30' Extended 
GOODS-South (aka the GEMS field, hereafter ECDFS; Rix 
et al. 2004; Castellano et al. 2010), the 30' x 30' Extended 
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Figure 1. Comparison of SEDS 3.6 /xm depth and total area (solid circle) to 
other major Spitzer / IRAC extragalactic surveys (open circles). Points indicate 
either achieved or SENS-PET la point-source sensitivities for GOODS (Great 
Observatories Origins Deep Survey), the EGS (Extended Groth Strip), E-CDFS 
(Extended Chandra Deep Field South), SpUDS (Spitzer Public Legacy Survey 
of UKIDSS Ultra-Deep Survey), SCOSMOS (Spitzer Deep Survey of HST 
COSMOS 2-Degree ACS Field), SERVS (Spitzer Extragalactic Representative 
Volume Survey), BCS (Blanco Cluster Survey), SWIRE (Spitzer Wide-area 
Infrared Extragalactic Survey), the FLS (Spitzer First-Look Survey), SDWFS 
(Spitzer Deep, Wide-Field Survey), the SSDF (SPT-Spitzer Deep Field), 
S-CANDELS (Spitzer-CANDELS), the UDF (Ultra-Deep Field), SIMPLE 
(the Spitzer IRAC/MUSYC Public Legacy in E-CDFS), SpIES (Spitzer-IRAC 
Equatorial Survey), and SPLASH (Spitzer Large- Area Survey with Hyper- 
Suprime-Cam). 
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GOODS-North (HDFN; Giavalisco et al. 2004; Wang et al. 
2010; Hathi et al. 2012; Lin et al. 2012), the 50' x 50' UKIDSS 
Ultra-Deep Survey (UDS; aka the Subaru/XMM Deep Field, 
Ouchi et al. 2001; Lawrence et al. 2007), a 10' x 60' region 
within the Extended Groth Strip (EGS; Davis et al. 2007; Bielby 
et al. 2012), and a 10' x 60' strip within the Ultra Vista deep 
survey of the larger COSMOS field (Scoville et al. 2007b; 
Koekemoer et al. 2007; McCracken et al. 2012). 

The five fields are distributed in ecliptic longitude and decli- 
nation to permit ground-based followup from both hemispheres. 
Each field has its own combination of strengths; no two are 
alike but each has an extensive suite of supporting data from 
both ground- and space-based observatories throughout the op- 
tical and near-IR bands. The deep near-IR imaging taken with 
Hubble Space Telescope ( HST) / WFC3 by Cosmic Assembly 
Near-Infrared Deep Extragalactic Legacy Survey (CANDELS; 
Grogin et al. 2011; Koekemoer et al. 2011) are an especially 
important component of these supporting data, because they 
cover a portion of each SEDS field in multiple infrared bands at 
significantly higher spatial resolution than IRAC. Some fields 
also have additional IRAC data that were not used either because 
they were relatively shallow or covered only limited areas. Notes 
on the IRAC observations of individual SEDS fields follow (all 
IRAC data sets used in this paper are listed in Table 1): 

The first IRAC observations of the UDS field were the 
shallow coverage obtained by SWIRE (4 x 30 s exposures, 
reaching 23.7 mag (lcr) at 4.5 /xm; Lonsdale et al. 2003, 2004; 
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Table 1 

The Five SEDS Fields 

Field PID a 

Epoch 3.6 fim. 

, 4.5 /cm BCDs Used 

UDS (2:18:00, -5:10:17, 0.32 deg 2 ) 

142 

2004 Jul 27-28 

548, 597 b 

40021 

2008 Jan 26-29 

3640, 3457 

61041 

2009 Sep 8-23 

5255, 5256 

61041 

2010 Feb 13-Mar 2 

5328, 5328 

61041 

2010 Sep 22-Oct 13 

5436, 5436 

ECDFS (3:32:20, -27:37:20, 0.35 deg 2 ) 

81 

2004 Feb 16 

167, 146 

194 

2004 Feb 8-16 

1724, 1723° 

194 

2004 Aug 12-18 

1632, 1632 c 

20708 

2005 Aug 19-23 

1943, 1872 

20708 

2006 Feb 6-1 1 

1899, 1944 

30866 

2007 Feb 15 

1200, 1080 

60022 

2010 Sep 20-Oct4 

4752, 4588 

60022 

201 1 Mar 26-Apr 3 

4596, 4752 

60022 

2011 Oct 10-Oct 20 

4717, 4552 

70204 

201 1 Mar 17-Apr 7 

5184,5128 

COSMOS (10:00:30, +2:10:00, 0.19 deg 2 ) 

20070 

2005 Dec 30-2006 Jan 2 

1259, 1253 

61043 

2010 Jan 25-Feb 4 

3672, 3672 

61043 

2010 Jun 10-28 

3164, 3140 

61043 

201 1 Jan 30-Feb 6 

3180, 3196 

HDFN (12:36:12, +62:14:12, 0.25 deg 2 ) 

81 

2004 May 26-27 

215, 178 

169 

2004 May 16-26 

2609, 2609 c 

169 

2004 Nov 17-25 

2447, 2447 c 

169 

2005 Nov 25 

114, 1 14 c 

20218 

2005 Nov 28-Dec 9 

200, 200 

20218 

2006 Jun 2-3 

200, 200 

61040 

2010 May 12-29 

4895, 4896 

61040 

201 1 Feb 28-Mar 13 

5440, 5440 

60140 

201 1 May 22-Jun 2 

5208, 4896 

EGS (14:19:38, +52:25:47, 0.35 deg 2 ) 

8 

2003 Dec 21-28 

988, 969 c 

8 

2004 Jun 28-Jul 3 

1027, 989° 

8 

2006 Mar 28-29 

117, 24 c 

41023 

2008 Jan 24-25 

726, 726 

41023 

2008 Jul 21-23 

726, 726 

61042 

2010 Feb 5-16 

4056, 4056 

61042 

2010 Aug 4-19 

4021,4056 

61042 

201 1 Feb 10-22 

3970, 4048 


Notes. SEDS field positions and areas. Areas are those covered in both IRAC 
channels to a depth of at least 10 ks by the projects listed. 
a Spitzer Program Identification Number. 
b 30 s frames. 
c 200 s frames. 

PID 142). Later, to match the depths of the newly-available 
near-IR photometry from the UKIRT Infrared Deep Sky Survey 
(UK1DSS; Lawrence et al. 2007) this field was surveyed a 
second time with IRAC as part of the Spitzer Public Legacy 
survey of the UKIDSS (SpUDS; Caputi et al. 2011). SpUDS 
covered a roughly 1 deg 2 area in all four of the IRAC bands 
available during the cryogenic mission to significantly deeper 
levels (e.g., reaching 24.7 mag (ler) at 4.5 /urn). 

The SEDS observations of the UDS were carried out in 
three epochs from 2009 to 2010 (Table 1). A 0.40 deg 2 region 
within SpUDS was covered during each SEDS visit with a 
total exposure time of roughly 72 x 100 s in each pixel of 


the resulting mosaic (Figure 2). When the full-depth UDS 
mosaic was assembled, it incorporated all the coextensive 
SWIRE and SpUDS exposures to reach the survey goal of 
12 hr total integration per pixel. The full-depth SEDS mosaics 
and coverage maps of the UDS are shown in Figures 3-6. 
The SEDS IRAC coverage is irregular because a portion 
of the northeast corner was avoided to prevent the bright 
stars there from introducing strong short-term latents into the 
exposures. A small area to the south was instead observed (the 
peninsula visible in Figures 5 and 6) because it was known 
to contain z — 6 Lya emitting galaxies. Users should be 
aware of this selection bias when using SEDS data from the 
peninsula. 

The ECDFS has been observed by a number of Spitzer pro- 
grams starting at the very beginning of the mission because of the 
compelling depth and wavelength coverage of complementary 
observations, including HST / ACS coverage from GEMS (Rix 
et al. 2004). The IRAC observations have been taken at a variety 
of depth/area combinations. The earliest IRAC mosaics come 
from the Great Observatories Origins Deep Survey (GOODS) 
project, and cover a small region to roughly 26 mag (3er) at 3.6 
and 4.5 /im. In Spitzer Cycle 2, a much larger 1600 arcmin 2 
region surrounding the GOODS region was covered by 
the IRAC/MuSYC Public Legacy Survey in the Extended 
CDF-South (SIMPLE; Damen et al. 2011). An additional pro- 
gram (Labbe et al. 2012) was carried out to get deep coverage of 
the two UDF NICMOS parallel fields with IRAC, resulting in an 
additional ~27 hr integration but covering a very small region, 
just slightly larger than a single IRAC field of view (FOV). 

The SEDS imaging, which covers a subregion of SIMPLE, 
was carried out in three separate visits in 2010 and 2011. 
The SEDS coverage was designed to avoid the GOODS IRAC 
pointings, which were already deeper than SEDS in any case. 
Because of the pre-existing coverage, SEDS only obtained 
roughly 108 x 100 s exposures at each pointing, this amount 
being sufficient to reach the target 12 hr total integration time. 
Finally, a Cycle 7 program (PID 70204) was carried out to 
image a very small region of the ECDFS but very deeply, 
accumulating an additional 75 hr in dithered 100 s exposures 
in approximately two IRAC FOVs within SEDS. The full-depth 
SEDS mosaics and coverage maps, which are combined with 
data from GOODS, SIMPLE, and PID 72004, are shown in 
Figures 7-10. 

The COSMOS field subtends roughly 2 deg 2 overall, located 
so as to take advantage of a local minimum in the far-infrared 
background (0.9 MJy sr ' 1 at 100/iiri). High-resolution HST/ 
ACS F814W imaging is available for the entire field (Scoville 
et al. 2007a; Koekemoer et al. 2007). COSMOS was first 
imaged by Spitzer in all four IRAC bands by S-COSMOS 
(Sanders et al. 2007). S-COSMOS consisted of dithered/ 
overlapping 12 x 100 s exposures, achieving sensitivities of 
23.9 and 23.3 mag (5a) at 3.6 and 4.5 /im, respectively, over 
2.3 deg 2 . 

The SEDS observations in S-COSMOS were arranged to 
cover a narrow strip roughly 10' x 1° oriented N-S. The coverage 
was sited to coincide with very deep ground-based YJHK S 
coverage obtained by McCracken et al. (2012) from the VISTA 
telescope (UltraVISTA). The SEDS data were accumulated in 
three visits in 2010 and 2011 (Table 1). The full-depth SEDS 
COSMOS mosaics and coverage maps, which incorporate 
the S-COSMOS exposures, are shown in Figures 1 1-14. The 
deep coverage is bordered by a relatively wide “crust” of 
intermediate-depth coverage. The crust arises from the tight 
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Figure 2. Cumulative area coverage as a function of exposure time for SEDS, including other, earlier observations (Table 1). The solid and dotted lines correspond 
respectively to the 3.6 and 4.5 /im bands. Each panel shows the coverage for the first SEDS epoch (lower curves) together with that for the full-depth mosaics that 
combine all data (Table 1) for each field. Panel (f) shows the coverage summed over all five SEDS fields. The nominal SEDS depth was 12 hr (43.2 ks); source 
extraction and photometry was only performed on regions having at least 10 ks depth. 


spacecraft roll angle constraints that prevail for this low-ecliptic 
latitude field. 

The HDFN was first observed by IRAC in Spitzer Cycle 1 by 
the GOODS team (Giavalisco et al. 2004), who used a narrow 
and deep configuration to cover the multiband data originally 
obtained with HST (Williams et al. 1996). Those first IRAC 
observations consisted of dithered 200 s exposures arranged 
in a pair of 2 x 2 maps with a deeper region of overlapping 
coverage in the center. Subsequently, a wider 300 arcmin 2 region 
surrounding the GOODS coverage was observed in Cycle 2 with 
100 s frames, reaching 24.5 mag at 3.6 /rm. 

SEDS was the next program to observe this region, covering it 
three times in 2010 and 2011. Each SEDS visit covered a region 
roughly 30' x 30' overlapping with the flanking-field coverage 
from Cycle 2 but avoiding the GOODS observations. The 
full-depth SEDS mosaics and coverage maps (Figures 15-18) 
incorporate the data from both previous programs. The coverage 
matches the footprint of the JHK survey conducted with CFHT / 
WIRCAM (Wang et al. 2010; Lin et al. 2012). 

The EGS, like the other four regions surveyed by SEDS, 
benefits from a substantial array of multiwavelength imaging 
(Davis et al. 2007). The first substantial Spitzer/IRAC imaging, 
reaching 50% completeness at 23.4 p , Jy at 3.6 /im, was carried 
out by Barmby et al. (2008) along a narrow 2? 3 x 12' strip 
placed to coincide with preexisting imaging acquired by HST s 
Advanced Camera for Surveys (ACS). These observations 


(PID 8, Fazio, Rieke, & Wright) were obtained in two separate 
visits in 2003 and 2004. They were augmented by two additional 
IRAC visits in 2008 (PID 41023, PI: Nandra) that broadened 
the central portion (~1 deg) of the long strip so as to provide 
much better overlap between the IRAC coverage and very deep 
(800 ks) X-ray imaging by Chcindra. 

The SEDS observations of the EGS cover the central degree, 
coincident with the IRAC imaging from both Barmby et al. 
(2008) and the additional observations from Cycle 4 as well 
as the multiband HST imaging from CANDELS. The SEDS 
imaging was carried out in three visits during 2010 and 2011 
(Table 1). The full-depth SEDS mosaics and coverage maps 
incorporate all coextensive data from the earlier programs. They 
are shown in Figures 19-22. 

3. OBSERVATIONS AND DATA REDUCTION 
3.1. Mapping Strategy 

The SEDS observing strategy was governed by the need to 
maximize the area covered with deep integrations in each of 
the five fields. It was also necessary to cover each field fully on 
separate visits so as to enable the planned variability science. 
Each field was visited three times; the visits (referred to below as 
epochs) were separated by six months except that the intervals 
for HDFN were nine and three months. Table 1 lists the epochs 
for each field. For each field, the coverage was very similar in 
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Figure 3. The full-depth SEDS 3.6 (tm mosaic in the UDS field (three coadded 
epochs, plus the coextensive imaging taken during the cryogenic mission 
and listed in Table 1). The image stretch ranges from —0.002 (white) to 
0.05 MJy sr" 1 (black). 
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Figure 4. As Figure 3, but showing the SEDS 4.5 /am mosaic of the UDS. The 
stretch ranges from —0.002 to 0.05 MJy sr -1 . 
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(black); the deepest coverage in this field is 65.1 ks. The total area observed 
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(A color version of this figure is available in the online journal.) 
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Figure 7. As Figure 3, but showing the SEDS IRAC 3.6 (tm mosaic of the 
ECDFS. The stretch ranges from 0 to 0.05 MJy sr^ 1 . 
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Figure 8. As Figure 7, but showing the full-depth mosaic of the ECDFS 
including all SEDS and cryogenic imaging by IRAC at 4.5 fim. The stretch 
ranges from 0 to 0.05 MJy sr" 1 . 
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Figure 9. The total IRAC 3.6 /im coverage map in the ECDFS including all 
data from SEDS as well as the cryogenic mission. The stretch is logarithmic and 
ranges from 0 (white) to 300 ks. The total area observed to at least 10 ks depth is 
0.34 deg 2 . The area surveyed by the HST MCT project CANDELS with WFC3 
and ACS is indicated in red. The deep PID 70204 imaging appears as adjacent, 
very dark squares inside the CANDELS field. The deepest IRAC integration is 
well over 100 hr. 

(A color version of this figure is available in the online journal.) 
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Figure 10. As Figure 9, but showing the full-depth IRAC 4.5 /xm coverage map 
in the ECDFS. The stretch is logarithmic and ranges from 0 (white) to 300 ks. 
The total area observed to at least 10 ks depth is 0.34 deg 2 . 
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Figure 11. The total IRAC 3.6 /im mosaic in the COSMOS field including all 
observations from SEDS and the cryogenic mission (Table 1). The image stretch 
ranges from —0.002 (white) to 0.05 MJy sr _1 (black). 
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Figure 12. As Figure 1 1 , but showing the full-depth 4.5 gm mosaic in the 
COSMOS field. The stretch ranges from —0.002 to 0.05 MJy sr -1 . 


epochs 1 and 3, but the coverage for epoch 2 differed because 
of the offset placement of the two IRAC detectors within the 
Spitzer focal plane and the constraints imposed by the spacecraft 
roll angle. 

Each SEDS epoch accumulated 4 hr integration time per 
pointing or less when there was pre-existing coverage of a field. 
Each Astronomical Observation Request (AOR), representing 
a continuous sequence of Spitzer observations, looked at each 
pointing position exactly once and covered an entire field. 
For efficiency, each position in the field was observed for 2-4 
frames of 100 s each before moving the telescope to the next 
position. There was no dithering built into individual AORs, 
but AOR origins were offset by ~1' in a Rouleaux pattern. 
Thus each field required 36-72 observations at each pointing 
(minus pre-existing coverage), and these were spread over the 
nearly the entire range of each epoch. Therefore each field 
was sampled with near-uniform time coverage and point-spread 
function (PSF). Individual AORs had different spacings between 
pointings and pointings observed in different time orders. The 
AORs thus sampled each sky position at a variety of positions 
on the arrays and also ensured that those positions included 


33 In some cases, grouped AORs were needed instead of a single one, but the 
effect was still to cover the entire field in a time interval of less than a day. 


both near and far spacings. These characteristics minimize data 
artifacts. 

Two of the fields, EGS and COSMOS, are long and narrow 
(~1° x 10')- The map pattern for these was two IRAC FOVs 
(each 5') in width by a variable number of FOVs in length. Vis- 
ibility constraints prevented the COSMOS strip from aligning 
with the IRAC arrays, and the width mapped to full depth in 
both FOVs was ~8'. The other fields (UDS, HDFN, ECDFS) 
are approximately square, and the basic map grid was rectangu- 
lar but with varying spacings. A few positions in the centers of 
the HDFN and ECDFS fields already possessed deep coverage 
(>23 hr), and pointings that would duplicate this coverage (four 
in ECDFS, two in HDFN) were omitted. The UDS field has 
bright stars in its northeast corner, and coverage was designed 
to exclude these stars. 

3.2. Data Reduction 

To the maximum extent possible, identical reduction proce- 
dures were applied to all SEDS data so as to facilitate subse- 
quent variability studies and to ensure a uniform data quality 
throughout the SEDS project. The data reduction was based 
on version SI 8.8.0 of the IRAC Corrected Basic Calibrated 
Data (cBCD) exposures. All 3.6 and 4.5 /xm cBCD frames were 
object-masked and median-stacked on a per-AOR basis; the 
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Figure 13. The total SEDS 3.6 /cm coverage map of COSMOS corresponding 
to the mosaic shown in Figure 11. The linear stretch ranges from 0 (white) to 
70 ks (black). The deepest coverage is 58.4 ks. The total area observed with at 
least 10 ks depth is 0.22 deg 2 . The areas surveyed with the WFC3 and ACS 
instruments by CANDELS are outlined respectively in red and blue. 

(A color version of this figure is available in the online journal.) 


resulting stacked images (presumed to represent blank sky) 
were visually inspected and subtracted from individual cBCDs 
within each AOR to eliminate long-term residual images arising 
from prior observations of bright sources. The sky-subtracted 
cBCDs were then examined individually and processed using 
custom software routines to correct column-pulldown effects 
associated with bright sources. The code, known as the “Warm- 
Mission Column Pulldown Corrector,” is publicly available at 
the Spitzer Science Center. ' 4 

After these preliminaries, the SEDS exposures and the coin- 
cident IRAC imaging from earlier programs (Section 2) were 
mosaiced together into 3.6 and 4.5 /im mosaics using IRACproc 
(Schuster et al. 2006). IRACproc implements the standard IRAC 
reduction software (MOPEX). The standard MOPEX error 
propagation techniques, which are based on the IRAC uncer- 
tainty images, were employed by IRACproc, with one signifi- 
cant exception. IRACproc augments the capabilities of MOPEX 
by calculating the spatial derivative of each image and adjusting 
the clipping algorithm accordingly. Pixels where the derivative 

34 http://ssc.spitzer.caltech.edu/archanaly/contributed/browse.html 
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Figure 14. As Figure 13, but showing the full-depth 4.5 (tm coverage map of 
the COSMOS field. The linear stretch ranges from 0 (white) to 70 ks (black); 
the deepest coverage in this field is 55 ks. The total area observed with at least 
10 ks depth is 0.21 deg 2 . 


is low (in the field) are clipped more aggressively than are pixels 
where the spatial derivative is high (point sources). This avoids 
downward biasing of point source fluxes in the output mosaics 
that might otherwise occur because of the slightly undersampled 
IRAC PSF. The software was configured to automatically flag 
and reject cosmic ray hits based on pipeline-generated masks to- 
gether with the adjusted sigma-clipping algorithm for spatially 
coincident pixels. 

In order to take advantage of the subpixel shifts of our 
mapping strategy, and to facilitate registration of the IRAC 
mosaics to the coextensive CANDELS imaging (O'.'Ob pixels; 
Grogin et al. 2011; Koekemoer et al. 2011), the mosaics 
were resampled to 0'.'6 pixel -1 . Thus each SEDS mosaic 
pixel subtends approximately one-fourth the area of the native 
IRAC pixel. The final mosaics’ tangent-plane projections were 
likewise aligned to those used by the CANDELS mosaics. 

The 40 final mosaics and coverage maps (including three 
epochs and a total coadd in each of the two IRAC bands for 
each of the five fields) are all available from the SEDS team 
Web site. 35 


35 http://www.cfa.harvard.edu/SEDS 
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Figure 15. Total SEDS IRAC 3.6 /xm mosaic in the HDFN field, with the 
observations from the warm mission included. The stretch ranges from 0 to 
0.05 MJy sr -1 . 
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Figure 16. As Figure 15, but for the 4.5 /xm mosaic. 


4. SOURCE EXTRACTION AND PHOTOMETRY 
4.1. Source Identification 


Figure 17. Total SEDS IRAC 3.6 /xm coverage map in the HDFN field. The 
logarithmic stretch ranges from 0 (white) to 2000; the deepest coverage in this 
field is roughly 160 ks. The red rectangle indicates approximately the area 
covered by the CANDELS observations. The total area observed with at least 
10 ks is 0.30 deg 2 . 

(A color version of this figure is available in the online journal.) 
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Figure 18. As Figure 17, but for the SEDS 4.5 fim coverage. The total area 
observed with at least 10 ks is 0.28 deg 2 . 


A first test of source extraction methods was made in 
the SEDS EGS field with SExtractor (ver. 2.5.0; Bertin & 
Arnouts 1996), a standard tool for these purposes (e.g., Lonsdale 
et al. 2003; Ashby et al. 2009). The software was configured 
identically to that used by Barmby et al. (2008) to photometer 
sources in their (shallower) EGS mosaic. Despite the fact that 


the SEDS EGS mosaics incorporate a factor of four longer 
integration time per pixel, the resulting SExtractor catalogs 
improved only marginally upon those of Barmby et al. In effect, 
source confusion imposed a limit on our ability to identify faint 
sources that SExtractor could not overcome. This behavior has 
been remarked on before (e.g., Sanders et al. 2007). We therefore 
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Figure 19. Total SEDS IRAC 3.6 /rm mosaic in the EGS held. The stretch 
ranges from 0 to 0.05 MJy sr _1 . 
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Figure 20. Total SEDS IRAC 4.5 pm mosaic in the EGS held. The stretch 
ranges from 0 to 0.05 MJy sr -1 . 
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used StarFinder (version 1.6f; Diolaiti et al. 2000) to identify 
IRAC sources in the SEDS fields. StarFinder is optimized 
for identification of blended sources in crowded wide-field 
adaptive-optics observations. StarFinder works by repeatedly 
fitting and subtracting a PSF to sources it identifies in an image. 
It begins with the brightest object. After the brightest source has 
been identified, fitted, and subtracted, the next brightest object 
is treated in the resulting residual image. StarFinder iterates 
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Figure 21. Total SEDS IRAC 3.6 i±m coverage map in the EGS field. The linear 
stretch ranges from 0 (white) to 700; the deepest effective coverage shown is 
roughly 20.3 hr. The total area observed with at least 10 ks is 0.26 deg 2 . The 
areas surveyed with the WFC3 and ACS instruments by CANDELS are outlined 
respectively in red and blue. 

(A color version of this figure is available in the online journal.) 
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Figure 22. Total SEDS IRAC 3.6 pm coverage map in the EGS field. The linear 
stretch ranges from 0 (white) to 700; the deepest effective coverage shown is 
roughly 21.1 hr. The total area observed with at least 10 ks is 0.25 deg 2 . 
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this process, operating on progressively fainter sources, until a 
user-specified limiting sensitivity is reached. 

The SEDS catalogs were constructed in two steps. First, 
StarFinder was used to locate sources (even faint, blended 
ones). Second, a custom code was used to correct biases in 
the StarFinder photometry. 
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Table 2 

SEDS Aperture Corrections 


Field Aperture Diameters (arcseconds) 



2.4 

3.6 

4.8 

6.0 

7.2 

12.0 

3.6 n m 

UDS 

2.29 

1.52 

1.29 

1.20 

1.16 

1.08 

ECDFS 

2.16 

1.50 

1.28 

1.17 

1.14 

1.07 

COSMOS 

2.11 

1.45 

1.25 

1.16 

1.13 

1.06 

HDFN 

2.28 

1.51 

1.28 

1.19 

1.15 

1.07 

EGS 

2.26 

1.54 

1.30 

1.19 

1.16 

1.08 

4.5 (im 

UDS 

2.27 

1.57 

1.33 

1.21 

1.18 

1.10 

ECDFS 

2.19 

1.49 

1.27 

1.18 

1.15 

1.07 

COSMOS 

2.15 

1.50 

1.27 

1.17 

1.14 

1.07 

HDFN 

2.22 

1.52 

1.29 

1.18 

1.14 

1.08 

EGS 

2.24 

1.50 

1.28 

1.19 

1.16 

1.08 


Notes. Aperture corrections (multiplicative factors) derived from empirical PSFs 
constructed using at least 10 stars in each of the SEDS fields. The corrections 
are defined so that when applied, they convert SEDS aperture photometry to 
total 24" aperture diameter magnitudes, as used for Spitzer calibration stars. 
All photometry in the SEDS catalogs presented in this work has been aperture- 
corrected using these values. 


The source-location step used the full-depth SEDS mosaics 
masked to exclude regions of shallow coverage. This exclusion 
of relatively noisy low-coverage regions — virtually all located 
on the periphery of the SEDS footprints — was done to prevent 
them from artificially inflating the rms-threshold criteria needed 
to search deeply for faint objects in the areas of highest signal/ 
noise. The 3.6 and 4.5 p m mosaics were masked separately, i.e., 
the coverage was not initially required to be coextensive. 

The source-location process began with the creation of 
empirical PSFs for the SEDS fields. Isolated, unsaturated stars 
(between magnitudes 10 and 15, depending on the field and 
the coverage) were inspected (to exclude binaries) and then 
scaled to a common, fiducial intensity. Objects within 64 mosaic 
pixels (38'.'4) of each PSF star were masked, and the resulting 
scaled images were median stacked to suppress artifacts. The 
halos of the resulting composite, high-dynamic-range PSFs were 
then smoothed slightly at large radii to suppress noise. They 
were also masked at radii beyond 64 pixels. Because the scale 
factors applied prior to the median stacking can be large for 
relatively faint stars, the scaling can lead to spurious negative- 
valued features that create obvious artifacts in the residual 
images described below. We therefore masked negative-valued 
pixels in the composite PSFs to zero. All the PSFs were then 
normalized to unity integrated flux, per standard StarFinder 
protocol. Correction factors were then estimated by measuring 
the flux from each PSF star within a series of apertures covering 
a range of diameters (Table 2). 

With the cleaned and normalized PSFs, iterative source 
identification was carried out as described above on the masked 
mosaics. To prevent bright sources from distorting the estimates 
of backgrounds nearby, StarFinder was configured to estimate 
the backgrounds with a grid spacing equal to 72 times the 
FWHM of the PSF. It was also set to estimate backgrounds under 
detected sources. A spacing constraint was imposed, requiring 
that sources be separated by at least the half of the PSF width 
(defined as the half-width at half- maximum or HWHM; 9 ) 
in order to be counted as separate detections; blended sources 
lying closer to each other were counted as a single object. 


StarFinder computes the rms of the source-free pixels to estimate 
the significance of source detections. In the first StarFinder 
run, sources were counted as detections if they exceeded 5 a 
significance relative to the estimated rms variations. After the 
first run had removed the 5a sources, StarFinder was run a 
second time but with a 3a threshold on the residual images, i.e., 
a new (lower) rms was computed and then applied to science 
mosaics from which all previously-detected sources had been 
removed. Finally, StarFinder was run a third time on the source- 
subtracted residual from the second iteration, again with a 3a 
detection criterion. Figure 23 shows StarFinder residual images 
for a typical small subregion (in this case, within the EGS). 
Although bright and extended sources have significant residuals, 
those for the more abundant faint sources are sufficiently small 
that the residual image is far from being confused in either 
band. The SEDS catalogs consist of the source positions and 
magnitudes determined in all three StarFinder runs, modulo the 
photometric corrections described below. 

4.1.1. Photometric Corrections 

StarFinder is very efficient for identifying sources, but it 
measures only PSF-fitted photometry. While this is perfectly 
acceptable for some sources (Milky Way stars, QSOs), it is 
inadequate for galaxies because they are seldom point sources. It 
was therefore necessary to convert the PSF-fitted photometry to 
aperture photometry in order to measure total fluxes for galaxies. 

To generate the SEDS aperture photometry we used a custom 
code that operated on the StarFinder catalogs and residual 
images. The code first re-created the original appearance of 
each source in turn by re-inserting the scaled/fitted PSF at the 
position of that object, effectively adding in the fitted flux to 
any residuals that remained. Thus each source was re-created 
at its original position with its original flux, with the PSF fits 
to all other sources removed. The code then re-measured the 
flux of this source in a suite of apertures (Table 2). Backgrounds 
were estimated within annuli centered on the sources but outside 
the photometric apertures. The procedure was repeated for all 
detected sources. This approach had two significant advantages. 
First, it greatly improved the background estimates because to a 
good approximation all neighboring sources had been removed. 
Second, it accounted for the fact that galaxies are not generally 
point sources; aperture photometry requires no assumptions 
about the surface brightness distributions. 

The multiaperture photometry was corrected to total magni- 
tudes using the aperture corrections measured separately for the 
empirical PSFs in each band and field. All SEDS catalog entries 
include these aperture corrections, which are given in Table 2 
and Figure 24. 

4.2. Completeness and Depth Simulations 

SEDS survey completeness and reliability were assessed 
with the standard Monte Carlo approach, using photometry of 
artificial sources. To constrain the properties of the required 
artificial sources for SEDS , we used a O'.' 5 search radius to match 
the SEDS COSMOS catalog to the coextensive F160W catalog 
made available by the CANDELS team (Grogin et al. 2011; 
Koekemoer et al. 2011). Figure 25 shows the measured widths 
(major axes only) of the CANDELS counterparts in the F160W 
filter. The great majority of SEDS IRAC-detected objects are 
significantly smaller than the IRAC PSF; they are in effect point 
sources at 3.6 /zm. At the magnitudes of interest to SEDS, i.e., 
objects fainter than 23 AB mag, virtually all the IRAC detections 
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Figure 23. StarFinder fits and subtracts sources from the crowded SEDS mosaics (left) to create residual images (right) that are relatively free of contaminating flux, 
making it possible to identify and photometer faint IRAC sources. The linear stretch ranges from —0.01 MJy sr -1 (white) to 0.04 MJy sr 1 (black) in all four panels in 
order to show the details of faint sources. Upper left panel: A 5' x 5' region in the 3.6 (tm EGS mosaic. Upper right panel: The same field after all StarFinder-detected 
sources have been scaled and subtracted. Lower left panel: The same field as it appears in the 4.5 (tm mosaic. Lower right panel: The 4.5 /am residual image. 


are point sources. We therefore used only point sources in our 
completeness and depth simulations. 

A minimum of 10 5 artificial point sources were introduced 
into each SEDS mosaic. The number of artificial sources 
inserted simultaneously, however, was kept small (less than 2% 
of the number of detected IRAC sources) so as not to artificially 
aggravate source confusion effects. Typically, only 50-250 
artificial sources were introduced at a time. Approximately the 
same number of sources was inserted in each 0.5 mag interval 
between 18 and 27 mag. They were placed at random locations 
throughout the portions of the science mosaics having at least 
10 ks integration time to sample all variations in large-scale 
structure present. The modified mosaics were then photometered 
with StarFinder in the same way the original detections had been 
made. 

We followed Barmby et al. (2008) and counted an artificial 
source as a valid detection if its measured flux was within 50% 
of its "true” flux and its measured position was within 1" of 
its a priori known position. Completeness was thus defined 
as the fraction of artificial sources over a range of apparent 
brightnesses in bins of width 0.5 mag extending from 100% 
completeness (at ~18 mag) down to zero at 26.5-27.0 mag. 
The uncertainties in the completeness estimates were inferred 


from Poisson statistics in each bin. As Figure 26 and Table 3 
show, the 3.6 /rm SEDS catalog is 80% and 50% complete at 
roughly 23.5 and 24.75 mag, respectively. The corresponding 
4.5 /im limits are 23.8 and 24.8 mag. Thus at magnitudes where 
the incompleteness correction becomes large (Figure 25), the 
choice of point sources for the simulations is appropriate. Some 
field-to-field variation is apparent in both bands. 

Roughly 70% of the faintest 3.6 /mi StarFinder detections 
(25.5-26.0 mag) in the COSMOS field have no F160W coun- 
terparts. The CANDELS F160W observations reached approx- 
imately 26.5 mag (3er), so genuine SEDS sources with typical 
colors ought to have been detected by CANDELS. The com- 
parison suggests that StarFinder generated spurious detections 
at faint levels. We therefore imposed a requirement that sources 
must be detected at both 3.6 and 4.5 /rm to be included in the 
SEDS catalogs. Formally, the requirement imposed was that the 
3.6)4111 detections match a 4.5 /im source within O'.' 5. Nearly 
100% of the sources brighter than 24.5 mag at 3.6 /im were 
detected in both SEDS bands. At levels fainter than 24.5 mag, 
however, the match fraction decreases with magnitude down to 
the SEDS survey limit. For this reason, Figure 26 and Table 3 
report the separate single-band completeness measurements for 
sources brighter than 24.5 mag but report (identical) two-band 
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Aperture Radius (arcsec) 


Figure 24. Aperture corrections (Table 2 factor minus one) in the five SEDS fields. Lines show empirical corrections derived from aperture photometry of the stacked 
PSFs of bright, unsaturated stars. Symbols represent the aperture corrections tabulated in the Spitzer Observers Manual for theoretical IRAC PSFs. 



Figure 25. The angular sizes of IRAC-detected sources in the COSMOS field. Points indicate the measured widths (major axes only) of CANDELS F160W counterparts 
to the SEDS 3.6 fim detections. The mean F160W FWHM and 1 <t uncertainties within bins of width 0.5 mag are indicated in red. Only a minority of such sources are 
extended compared to the FWHM of the IRAC 3.6 fim PSF (indicated by the blue dotted line). All sources fainter than 23.5 AB mag are effectively point sources for 
SEDS. 

(A color version of this figure is available in the online journal.) 


completeness estimates for all SEDS sources fainter than this 
level. 

The common practice (e.g., SExtractor) of basing photometric 
uncertainties on noise estimates made in source-free mosaic 
pixels can be problematic in some circumstances. The issue 


of most concern to SEDS is that the supposedly source-free 
pixels may well fall on faint sources lying just below the 
detection threshold; their fluxes may then artificially boost 
the uncertainty estimates. Also, sub-pixel sampling like that 
used for SEDS introduces correlated noise into the mosaics. 


13 


The Astrophysical Journal, 769:80 (27pp), 2013 May 20 


Ashby et al. 



AB Magnitude 

Figure 26. Completeness in the SEDS fields estimated by Monte Carlo simulations as described in Section 4.2. Solid triangles and squares correspond to the IRAC 
3.6 and 4.5 /im bands, respectively. The values are averages over the field areas having at least 10 ks exposure time, but completeness in any particular small region 
depends on the actual integration time achieved there. 


To avoid these complications the SEDS error estimates were 
based solely on photometry of the simulated sources with their 
known fluxes. We grouped the artificial sources into 0.5 mag bins 
and measured the offsets from the known fluxes as a function 
of magnitude. The standard deviations of the offsets were taken 
as the ler measurement uncertainties and the offsets themselves 
were taken as the measurement biases. The results are shown in 
Figure 27. 

The uncertainties are much more significant than the biases 
in all SEDS fields, at least for sources brighter than 26 mag. The 
uncertainty typically rises smoothly from roughly 0.03 mag for 
bright sources to about 0.25 mag at 26 AB mag. It combines 
contributions from three separate sources of uncertainty: (1) the 
3% uncertainty in the IRAC absolute calibration (Reach et al. 
2005), (2) photon noise, and (3) confusion noise. These are 
discussed further in Section 6. 1 . The empirical uncertainties are 
given in Table 4. 

The measurement bias is relatively small for sources brighter 
than 20 mag. It starts to grow rapidly at roughly the level 
where SEDS becomes 50% incomplete and becomes greater 
than 0.1 mag at 26 AB mag. This is consistent with a picture 
in which faint sources are increasingly difficult to deblend from 
on-average brighter neighbors. The contamination of the photo- 
metric apertures by these neighbors affects the photometry, even 
though the measurements were made in residual images. The 
measurement bias is given in Table 5. The photometry reported 
in the SEDS catalogs has been corrected to remove this bias. 


4.3. Verification 

Because SEDS extends to extremely faint levels (by current 
standards) and must cope with pervasive source confusion, we 
verified the SEDS measurements by comparing them to previ- 
ously published catalogs. We estimated the astrometric uncer- 
tainties by comparing to the USNOB and 2MASS Point Source 
Catalog. We verified our photometric calibration by comparing 
to coextensive surveys including S-COSMOS (Sanders et al. 
2007), SpUDS (version DR2), EGS (Barmby et al. 2008), and 
SIMPLE (Damen et al. 201 1). 

4.3.1. Astrometric Reliability 

To estimate the accuracy of SEDS astrometry, we compared 
the SEDS positions of bright but unsaturated sources to those 
in the Two Micron All Sky Survey (2MASS; Skrutskie et al. 
2006) Point Source Catalog. We performed a search within 
1" of the positions of IRAC sources to identify their 2MASS 
counterparts. The distributions of coordinate offsets for the 3.6 
and 4.5 jim sources are shown in Figure 28. The astrometric 
discrepancies are small compared to the size of a SEDS 
pixel: at 3.6 /rm, the mean difference (SEDS— 2MASS) was 
just —O'.' 002 ± O'.' 16 in right ascension and O'.' 004 ± O'.' 19 in 
declination. The corresponding mean offsets for the 4.5 jim 
sources are similar: — O'/ 005 ± 0'.' 1 6 in right ascension and 
O'.' 03 ± O'.' 1 8 in declination. The radial uncertainties are therefore 
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Table 3 

Completeness in the SEDS IRAC Catalogs 


AB Mag UDS ECDFS COSMOS HDFN EGS 

3.6 fim 


18.25 

1.00 

± 

0.03 

1.00 

± 

0.03 

1.00 

± 

0.05 

1.00 

it 

0.03 

1.00 

± 

0.03 

18.75 

1.00 

± 

0.03 

0.99 

± 

0.03 

1.00 

± 

0.05 

1.00 

± 

0.03 

0.99 

± 

0.03 

19.25 

0.99 

± 

0.02 

0.99 

± 

0.03 

0.99 

± 

0.02 

0.99 

± 

0.02 

0.99 

± 

0.03 

19.75 

0.99 

± 

0.02 

0.99 

± 

0.02 

0.99 

d= 

0.02 

0.99 

± 

0.01 

0.98 

± 

0.02 

20.25 

0.98 

± 

0.02 

0.98 

± 

0.02 

0.98 

± 

0.02 

0.98 

± 

0.01 

0.97 

± 

0.02 

20.75 

0.97 

± 

0.02 

0.97 

± 

0.02 

0.96 

± 

0.02 

0.96 

± 

0.01 

0.96 

± 

0.02 

21.25 

0.95 

± 

0.02 

0.96 

± 

0.03 

0.95 

± 

0.02 

0.95 

± 

0.01 

0.94 

± 

0.03 

21.75 

0.93 

± 

0.02 

0.94 

± 

0.03 

0.93 

± 

0.02 

0.92 

± 

0.01 

0.93 

± 

0.02 

22.25 

0.91 

± 

0.02 

0.92 

± 

0.03 

0.90 

± 

0.02 

0.90 

it 

0.02 

0.90 

± 

0.02 

22.75 

0.88 

± 

0.02 

0.89 

± 

0.02 

0.87 

± 

0.01 

0.86 

± 

0.02 

0.87 

± 

0.02 

23.25 

0.82 

± 

0.01 

0.85 

± 

0.01 

0.82 

± 

0.01 

0.81 

± 

0.01 

0.82 

± 

0.01 

23.75 

0.78 

± 

0.01 

0.79 

± 

0.01 

0.77 

± 

0.01 

0.75 

± 

0.01 

0.75 

± 

0.01 

24.25 

0.69 

± 

0.01 

0.73 

± 

0.01 

0.68 

± 

0.01 

0.67 

± 

0.01 

0.67 

± 

0.01 

4.5 fim 

18.25 

1.00 

± 

0.06 

0.99 

± 

0.06 

1.00 

± 

0.03 

1.00 

± 

0.07 

1.00 

± 

0.06 

18.75 

1.00 

± 

0.06 

1.00 

± 

0.06 

1.00 

± 

0.03 

0.99 

it 

0.07 

1.00 

± 

0.06 

19.25 

0.99 

± 

0.02 

1.00 

± 

0.03 

0.99 

± 

0.02 

0.99 

± 

0.03 

0.99 

± 

0.03 

19.75 

0.99 

± 

0.02 

0.99 

± 

0.02 

0.99 

± 

0.02 

0.99 

± 

0.02 

0.99 

± 

0.02 

20.25 

0.98 

± 

0.02 

0.99 

± 

0.02 

0.98 

it 

0.02 

0.99 

± 

0.02 

0.99 

± 

0.02 

20.75 

0.98 

it 

0.02 

0.98 

± 

0.02 

0.97 

± 

0.02 

0.97 

± 

0.02 

0.97 

± 

0.02 

21.25 

0.96 

± 

0.02 

0.96 

± 

0.02 

0.96 

± 

0.02 

0.96 

± 

0.02 

0.95 

± 

0.02 

21.75 

0.94 

± 

0.02 

0.95 

± 

0.02 

0.94 

± 

0.02 

0.95 

± 

0.02 

0.94 

± 

0.02 

22.25 

0.93 

± 

0.02 

0.93 

± 

0.02 

0.91 

± 

0.02 

0.93 

± 

0.02 

0.94 

± 

0.02 

22.75 

0.90 

± 

0.01 

0.91 

± 

0.01 

0.88 

± 

0.01 

0.91 

d= 

0.02 

0.90 

± 

0.01 

23.25 

0.86 

± 

0.01 

0.87 

± 

0.01 

0.84 

± 

0.01 

0.85 

± 

0.01 

0.86 

± 

0.01 

23.75 

0.82 

± 

0.01 

0.83 

± 

0.01 

0.79 

± 

0.01 

0.81 

± 

0.01 

0.81 

± 

0.01 

24.25 

0.74 

± 

0.01 

0.76 

± 

0.01 

0.70 

± 

0.01 

0.74 

± 

0.01 

0.75 

± 

0.01 

3.6 and 4.5 fim 

24.75 

0.46 

± 

0.02 

0.56 

± 

0.03 

0.45 

± 

0.02 

0.52 

± 

0.01 

0.47 

± 

0.02 

25.25 

0.19 

± 

0.01 

0.30 

± 

0.02 

0.21 

± 

0.01 

0.31 

± 

0.01 

0.20 

± 

0.01 

25.75 

0.03 

± 

0.006 

0.11 

± 

0.01 

0.04 

± 

0.005 

0.11 

± 

0.008 

0.05 

± 

0.01 

26.25 

0.003 

± 

0.002 

0.008 

± 

0.003 

0.007 

± 

0.002 

0.03 

± 

0.004 

0.007 

± 

0.003 


Notes. Completeness estimates for the SEDS fields. The magnitudes correspond to the centers of bins of width 0.5 mag in which the 
completeness was estimated. The completeness is unity at brighter magnitudes than those listed. At magnitudes ^24.5, completeness 
estimates require detection in both IRAC bands. 


0'.'24 (ler) relative to 2MASS, or about 0.4 of a SEDS mosaic 
pixel. 

We also compared the SEDS positions to their visible- 
wavelength counterparts in the USNOB 1 .0 catalog (Monet et al. 
2003). We matched the positions in each of the five SEDS fields 
to those in USNOBl.O using a 0'.'5 search radius. Only USNOB 
sources having positions known to within O'.' 25 were matched. 
The mean offsets between USNOB and SEDS were insignificant 
for the UDS, ECDFS, and COSMOS fields. Significant offsets 
were found for the HDFN and EGS fields. The offsets are 
quantified for each field in Table 6. Users of the SEDS catalogs 
should be aware of and account for these offsets when matching 
to the USNOB positions. 

4.3.2. SEDS Photometric Verification 

To verify the SEDS photometry, we matched the SEDS cata- 
logs to those of the four coextensive surveys mentioned above. 
In all cases, the matching was done using a 0'.'' 5 search radius, 
i.e., roughly twice the lcr astrometric uncertainty established 
in Section 4.3.1. Only SEDS sources brighter than the detec- 
tion limits of the respective surveys were used in the compari- 
son. For example, only sources brighter than the SpuDS faint- 
source detection thresholds (2.0 and 2.7 //Jy at 3.6 and 4.5 //m, 


respectively) were matched to the SpUDS catalog. Only SEDS 
sources brighter than 1 (1.7) /rJy at 3.6 (4.5) // m were matched 
to S-COSMOS. The limits for the comparisons in the EGS field 
were set to match the 50% completeness threshold of Barmby 
et al. (2008), i.e., 1.5 and 1.6 /rJy at 3.6 and 4.5 /rm, respec- 
tively. Only SEDS sources brighter than 24.64 and 24.22 mag 
(the SIMPLE 5 ct sensitivities at 3.6 and 4.5 \i m, respectively) 
were matched to the Damen et al. (2011) SIMPLE catalog. In 
addition, sources flagged as blended in the SIMPLE DR1 cata- 
log were excluded from the comparison. Because all SIMPLE 
sources brighter than 18 AB mag were flagged as blended, the 
comparison could only be carried out for objects fainter than 
this limit. The results are shown in Figure 29. 

The SEDS photometry is broadly consistent with the previ- 
ous measurements. The mean offset for faint sources varies from 
zero in some cases up to a maximum of ~0.07 mag for the SIM- 
PLE 4.5 //m band. The comparisons for bright SEDS sources 
show a somewhat different character. With the exception of the 
EGS 3.6 /rm band, SEDS sources brighter than 18 mag show a 
systematic offset that ranges up to a maximum of 0. 1 mag in the 
worst case (versus SpUDS; Figures 29(a) and (b)) in the sense 
that these sources are brighter in the SEDS catalog on average. 
This may be due to the different source types that predominate on 
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AB Magnitude 


Figure 27. SEDS measurement errors based on photometry of simulated sources as described in Section 4.2. The symbols indicate the measurement bias as a function 
of apparent magnitude. The solid and dashed lines indicate the 1 o measurement uncertainty at 3.6 and 4.5 fim, respectively. The lower limit of 0.03 mag on the 
measurement uncertainties reflects current knowledge of the uncertainty in the IRAC absolute calibration. 



Figure 28. Left panel: Astrometric offsets (SEDS— 2MASS) measured between unsaturated SEDS sources and their counterparts in the 2MASS Point Source Catalog. 
A 1" search radius was used to match roughly 2500 2MASS sources (~500 objects per SEDS field) with objects detected by IRAC at 3.6 /xm. A high density of 
matches in this space is indicated with black. Gaussian functions (red) were fitted to the one-dimensional histograms of the offset distributions to estimate the la 
uncertainties, which are O'.' 16 in right ascension and O'.' 19 in declination. Right panel: The distribution for matched 4.5 fim SEDS sources. The lcr uncertainties derived 
from fitting a Gaussian profile to the histograms of offsets are O'.' 16 in right ascension and O'.' 18 in declination. 

(A color version of this figure is available in the online journal.) 


opposite sides of the 18 mag threshold. Whereas sources fainter 
than 18 mag are typically galaxies, objects brighter than 18 mag 
are much more likely to be Galactic stars, i.e., point sources 
(Arendt et al. 1998; Fazio et al. 2004a; see also Figure 25). 


Applying identical aperture corrections to both populations most 
likely introduced a slight bias to the SEDS 18-21 mag photom- 
etry of a few percent. This bias will not be significant for faint 
sources, because they are pointlike (Figure 25). 
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Table 4 

Empirical Photometric Uncertainties for SEDS 


Mag 

UDS 

ECDFS 

COSMOS 

HDFN 

EGS 

3.6 fim 

16.25 

0.03 

0.03 

0.03 

0.03 

0.03 

16.75 

0.04 

0.03 

0.03 

0.04 

0.03 

17.25 

0.05 

0.05 

0.03 

0.04 

0.03 

17.75 

0.06 

0.03 

0.03 

0.04 

0.05 

18.25 

0.05 

0.04 

0.03 

0.07 

0.07 

18.75 

0.07 

0.06 

0.03 

0.06 

0.07 

19.25 

0.07 

0.04 

0.06 

0.05 

0.07 

19.75 

0.08 

0.07 

0.08 

0.08 

0.08 

20.25 

0.08 

0.08 

0.08 

0.08 

0.08 

20.75 

0.09 

0.09 

0.09 

0.09 

0.09 

21.25 

0.10 

0.09 

0.10 

0.10 

0.10 

21.75 

0.11 

0.10 

0.11 

0.11 

0.12 

22.25 

0.13 

0.11 

0.12 

0.12 

0.12 

22.75 

0.13 

0.13 

0.14 

0.14 

0.13 

23.25 

0.15 

0.14 

0.15 

0.16 

0.16 

23.75 

0.18 

0.17 

0.18 

0.19 

0.19 

24.25 

0.22 

0.21 

0.22 

0.22 

0.22 

24.75 

0.25 

0.24 

0.25 

0.25 

0.25 

25.25 

0.30 

0.28 

0.30 

0.30 

0.30 

25.75 

0.34 

0.33 

0.33 

0.34 

0.34 

4.5 /rm 

17.75 

0.03 

0.03 

0.03 

0.03 

0.03 

18.25 

0.03 

0.03 

0.04 

0.03 

0.03 

18.75 

0.03 

0.03 

0.06 

0.03 

0.03 

19.25 

0.07 

0.03 

0.07 

0.03 

0.03 

19.75 

0.07 

0.07 

0.08 

0.08 

0.07 

20.25 

0.07 

0.08 

0.08 

0.08 

0.08 

20.75 

0.08 

0.08 

0.09 

0.09 

0.09 

21.25 

0.09 

0.09 

0.10 

0.09 

0.09 

21.75 

0.10 

0.10 

0.11 

0.10 

0.10 

22.25 

0.11 

0.10 

0.11 

0.11 

0.11 

22.75 

0.13 

0.12 

0.13 

0.13 

0.12 

23.25 

0.14 

0.14 

0.15 

0.14 

0.14 

23.75 

0.18 

0.17 

0.18 

0.17 

0.17 

24.25 

0.21 

0.20 

0.22 

0.21 

0.21 

24.75 

0.25 

0.23 

0.25 

0.24 

0.24 

25.25 

0.30 

0.28 

0.30 

0.29 

0.30 

25.75 

0.34 

0.33 

0.34 

0.33 

0.34 


Notes. Empirically determined la photometric uncertainties (magnitudes) 
determined using the Monte Carlo simulations described in Section 4.2. An 
estimated 3% systematic error in the IRAC flux calibration is included and 
limits the uncertainties for bright sources. 

Several factors could account for the offsets. Changes in the 
IRAC calibration during the Spitzer mission could introduce 
changes no larger than 2%. The discrepancy already noted for 
the SIMPLE 4.5 p m band may in part be a result of the slightly 
different flux calibrations applied to the two SIMPLE epochs. 
However, a more likely cause for the discrepancies is the very 
different methods used to obtain the photometry in the first place. 
Whereas SEDS is based on a modified StarFinder approach, the 
SpUDS, S-COSMOS, EGS, and SIMPLE measurements are 
all based on SExtractor. In addition, the different surveys mea- 
sured their photometry in different apertures; the comparisons 
presented in Figure 29 used aperture diameters of 3" 8 (SpUDS), 
2" 8 (S-COSMOS), 4" 2 (EGS), and 4"0 (SIMPLE). Even small 
differences in applied aperture corrections could potentially in- 
troduce the offsets seen here. 

In summary, the comparison to earlier, coextensive sur- 
veys shows that SEDS faint-source photometry is generally 


Table 5 


Photometric Bias in SEDS Catalogs 

Mag 

UDS 

ECDFS 

COSMOS 

HDFN 

EGS 

3.6 fim 

16.25 

0.00 

0.00 

0.00 

0.00 

0.00 

16.75 

0.00 

0.00 

0.00 

0.00 

-0.01 

17.25 

0.00 

0.01 

0.00 

0.00 

-0.01 

17.75 

0.00 

0.00 

0.00 

0.00 

-0.01 

18.25 

0.00 

0.00 

0.00 

0.00 

-0.01 

18.75 

0.00 

0.00 

0.00 

0.01 

0.00 

19.25 

0.01 

0.00 

0.00 

0.01 

0.00 

19.75 

0.01 

0.01 

0.01 

0.01 

0.01 

20.25 

0.01 

0.01 

0.01 

0.01 

0.01 

20.75 

0.01 

0.01 

0.01 

0.02 

0.01 

21.25 

0.02 

0.01 

0.01 

0.02 

0.02 

21.75 

0.02 

0.02 

0.02 

0.02 

0.02 

22.25 

0.03 

0.01 

0.01 

0.02 

0.02 

22.75 

0.03 

0.01 

0.02 

0.03 

0.03 

23.25 

0.04 

0.02 

0.02 

0.03 

0.03 

23.75 

0.05 

0.02 

0.02 

0.05 

0.04 

24.25 

0.07 

0.04 

0.03 

0.06 

0.06 

24.75 

0.09 

0.03 

0.04 

0.08 

0.06 

25.25 

0.14 

0.07 

0.08 

0.12 

0.10 

25.75 

0.24 

0.13 

0.21 

0.19 

0.20 

4.5 fim 

16.25 

0.00 

0.00 

0.00 

0.00 

0.00 

16.75 

0.00 

0.00 

0.00 

0.00 

-0.01 

17.25 

0.00 

0.00 

0.00 

0.00 

-0.01 

17.75 

0.00 

0.00 

0.00 

0.00 

-0.01 

18.25 

0.00 

0.00 

0.00 

0.00 

-0.01 

18.75 

0.00 

0.00 

0.00 

0.00 

-0.01 

19.25 

0.01 

0.01 

0.01 

0.01 

-0.01 

19.75 

0.01 

0.01 

0.01 

0.01 

0.01 

20.25 

0.01 

0.01 

0.01 

0.01 

0.01 

20.75 

0.01 

0.01 

0.01 

0.01 

0.01 

21.25 

0.01 

0.01 

0.01 

0.01 

0.01 

21.75 

0.02 

0.01 

0.01 

0.02 

0.02 

22.25 

0.02 

0.01 

0.02 

0.02 

0.02 

22.75 

0.02 

0.02 

0.02 

0.02 

0.02 

23.25 

0.03 

0.02 

0.02 

0.02 

0.02 

23.75 

0.03 

0.02 

0.03 

0.03 

0.03 

24.25 

0.05 

0.03 

0.04 

0.04 

0.04 

24.75 

0.07 

0.04 

0.05 

0.05 

0.05 

25.25 

0.15 

0.08 

0.12 

0.08 

0.11 

25.75 

0.30 

0.21 

0.29 

0.19 

0.25 

Notes. Mean photometric bias 

in SEDS fields (magnitudes), 

determined 

empirically by the Monte Carlo simulations described in Section 4.2. The sense 
of the bias is that artificial sources are measured to be brighter, on average, than 

they were a 

priori known to be by the stated amounts. These biases have already 

been corrected in the catalogs presented here. 






Table 6 





Astrometric Offsets (SEDS-USNOB) 



Field 


AR.A. 

ADecl. 


Total 



(arcsec) 

(arcsec) 


(arcsec) 

UDS 


0.09 ±0.19 

0.00 ±0.18 


0.24 ±0.12 

ECDFS 


0.08 ± 0.20 

-0.03 ±0.19 


0.27 ±0.12 

COSMOS 


-0.02 ±0.19 

-0.11 ±0.17 


0.25 ±0.12 

HDFN 


0.18 ±0.24 

-0.28 ± 0.26 


0.44 ±0.19 

EGS 


-0.04 ± 0.25 

-0.27 ± 0.25 


0.39 ± 0.20 


Notes. Mean coordinate offsets measured between SEDS sources and their 
visible- wavelength counterparts in the USNOBl.O catalog (Monet et al. 2003). 
Total offsets refer to the mean absolute offsets between SEDS and USNOB 
positions. The stated uncertainties are the standard deviations of the offset 
distributions for matched sources. 
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Table 7 

SEDS Full-depth Source Catalog for the UDS Field 


Object 

R.A., Decl. 
(J2000) 


3.6 /xm Magnitudes 3 
4.5 /xm Magnitudes 


3.6 /xm Unc. b 
4.5 /xm Unc. 

3.6 (im Cvg c 
4.5 fi m Cvg e 

3.6 /xm Flag d 
4.5 /xm Flag f 

SEDS J021641. 84-051639.5 

34.17434-5.27763 

13.95 

13.94 

13.77 

13.68 

13.63 13.61 

13.57 

0.03 

149 

i 



14.11 

14.07 

13.97 

13.90 

13.88 13.87 

13.84 

0.03 

0 

3 

SEDS J021717. 25— 050411.5 

34.32189-5.06986 

13.98 

13.94 

13.66 

13.50 

13.40 13.36 

13.27 

0.03 

176 

1 



13.77 

13.69 

13.52 

13.42 

13.37 13.36 

13.32 

0.03 

152 

1 

SEDS J021737. 30— 045221.4 

34.40543 -4.87262 

14.00 

13.99 

13.87 

13.81 

13.78 13.76 

13.73 

0.03 

0 

3 



14.73 

14.64 

14.46 

14.35 

14.30 14.28 

14.24 

0.03 

0 

3 

SEDS J021756. 88— 050757.4 

34.48698-5.13260 

14.03 

14.02 

13.84 

13.76 

13.71 13.69 

13.65 

0.03 

187 

1 



14.08 

14.05 

13.96 

13.89 

13.86 13.86 

13.84 

0.03 

153 

1 

SEDS J021 81 8. 17—050009.9 

34.57571 -5.00276 

14.03 

13.97 

13.65 

13.45 

13.33 13.28 

13.17 

0.03 

145 

1 



13.75 

13.65 

13.43 

13.31 

13.25 13.23 

13.18 

0.03 

107 

1 

SEDS J021640. 10— 051740.5 

34.16709-5.29459 

14.11 

14.10 

14.00 

13.95 

13.92 13.91 

13.89 

0.03 

135 

1 



14.60 

14.57 

14.47 

14.40 

14.37 14.36 

14.34 

0.03 

0 

3 

SEDS J021740.09— 050249.6 

34.41705 -5.04712 

14.12 

14.11 

14.00 

13.94 

13.92 13.90 

13.88 

0.03 

237 

1 



14.20 

14.19 

14.13 

14.10 

14.08 14.08 

14.07 

0.03 

179 

1 

SEDS J021639. 71— 051836.8 

34.16547-5.31023 

14.20 

14.20 

14.13 

14.10 

14.08 14.08 

14.07 

0.03 

124 

1 



14.60 

14.58 

14.51 

14.46 

14.45 14.44 

14.43 

0.03 

0 

3 

SEDS J021725. 24— 051804.7 

34.35516-5.30130 

14.22 

14.21 

14.14 

14.10 

14.08 14.07 

14.05 

0.03 

226 

1 



14.34 

14.34 

14.32 

14.30 

14.29 14.29 

14.28 

0.03 

344 

1 

SEDS J021657. 24— 050801.5 

34.23849-5.13374 

14.22 

14.22 

14.14 

14.10 

14.08 14.08 

14.06 

0.03 

274 

1 



14.45 

14.44 

14.41 

14.39 

14.38 14.37 

14.37 

0.03 

210 

1 


Notes. The SEDS catalog of sources in the UDS field selected at both 3.6 and 4.5 gm as described in the text. 

a The PSF-fitted magnitude is given first, and the magnitudes given after are measured in apertures of 2'.'4, 3'.'6, 4'.'8, 6'.'0, 7'.'2, and 12'.'0 diameter. 
b Stated errors are la uncertainties for the 2'.'4 diameter aperture magnitudes. 

c Depth of coverage expressed in terms of the number of unsaturated 100 s 3.6 // ni IRAC frames that observed the source. 
d 3.6 //. m photometric quality flag, as described in Section 5. If the flag is nonzero, the photometry may be corrupted. 
e Depth of coverage expressed in terms of the number of unsaturated 100 s 4.5 fim IRAC frames that observed the source. 
f 4.5 //in photometric quality flag, as described in Section 5. If the flag is nonzero, the photometry may be corrupted. 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.) 


consistent with previous results even though the SEDS mea- 
surements were made using different software and methods and 
in significantly more crowded mosaics. Given the use of differ- 
ent photometric software with different parameterizations and 
the different depths of the data sets being compared, we regard 
the agreement between SEDS and the preexisting surveys as 
excellent. 

5. SEDS CATALOGS 

The SEDS IRAC catalogs are presented in Tables 7 
through 1 1 . In addition to the PSF-fitted magnitudes on which 
the source detection is based, they also list the IRAC positions 
with the 3.6 and 4.5 jrm photometry in six apertures of diame- 
ters 2"4, 3"6, 4"8, 6", 7'.'2, and 12". All the photometry has been 
aperture-corrected to total magnitudes assuming the sources are 
pointlike. It has also been adjusted to account for the empir- 
ically determined biases given in Table 5. In other words, all 
catalog magnitudes represent total magnitudes if the source is 
unresolved by IRAC, but the magnitudes of extended sources 
represent only the flux inside the given aperture. As shown in 
Figure 25, this will affect mainly galaxies with [3.6] < 21 mag. 
The catalogs also provide ler uncertainty estimates, but only for 
the 2" 4 diameter aperture photometry because time constraints 
made it impractical to simulate photometry for all the aperture 
diameters. Users are encouraged to use the 2" 4 aperture for 
photometry of faint sources. The other apertures are provided 
so that users can construct for themselves the curve of growth 
for large, extended sources. 

Users should be aware of some limitations of the SEDS source 
catalogs, which are described here. 


1. First and foremost, we do not present the multi-epoch 
photometry here. Catalogs corresponding to each visit 
to each SEDS field will be presented in a future 
contribution. 

2. Because of the large uncertainties (random and system- 
atic) associated with sources fainter than 26 AB mag, the 
SEDS catalogs have been cut at this magnitude. Specifi- 
cally, only sources detected at both 3.6 and 4.5 /im, and 
brighter than 26 mag in at least one band, appear in the cat- 
alog. Although SEDS is sufficiently deep to detect fainter 
sources, the signal-to-noise ratio falls below 3cr signifi- 
cance at 26 mag. Such faint sources are also subject to 
large systematic effects (e.g., the large systematic errors 
in the faint-source completeness estimates given in Ta- 
ble 3 dominate the uncertainties in the counts), and they 
are not reliably recovered in coextensive CANDELS HST / 
WFC3 F160W observations. Nonetheless, at 26 AB mag 
the SEDS catalogs reach fainter than any other catalogs 
of which we are aware for a survey of this size at these 
wavelengths. 

3. Confusion is a major issue for the SEDS catalogs, given 
the high source area density relative to the effective IRAC 
beam sizes at 3.6 and 4.5 /xm. We require SEDS sources 
to be farther apart than one IRAC HWHM (roughly 
(y/9), because sources closer to each other than this could 
not be reliably deblended and photometered. Thus close 
pairs or multiples may be cataloged by SEDS as single 
sources. 

4. Because the IRAC arrays view different areas of sky 
and their orientation is constrained by spacecraft visibility 
limits, the areas covered by the 3.6 and 4.5 jim detectors do 
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Table 8 

SEDS Full-depth Source Catalog for the ECDFS Field 


Object 

R.A., Decl. 
(J2000) 


3.6 fim Magnitudes 3 
4.5 fim Magnitudes 


3.6 fim Unc. b 
4.5 fim Unc. 

3.6 Atm Cvg c 
4.5 /Ltm Cvg e 

3.6 fim Flag d 
4.5 fim Flag f 

SEDS J033255. 60— 273039.4 

53.23165 -27.51095 

14.43 

14.38 

14.27 

14.21 

14.19 14.18 

14.15 

0.03 

0 

3 



14.88 

14.82 

14.70 

14.62 

14.58 14.57 

14.54 

0.03 

0 

3 

SEDS J033338. 65— 275033.7 

53.41103 -27.84270 

14.48 

14.44 

14.31 

14.25 

14.22 14.21 

14.18 

0.03 

0 

3 



15.05 

15.00 

14.85 

14.75 

14.71 14.70 

14.67 

0.03 

130 

0 

SEDS J033232. 62— 273754.0 

53.13592-27.63166 

14.49 

14.41 

14.27 

14.23 

14.21 14.20 

14.15 

0.03 

0 

3 



14.86 

14.81 

14.63 

14.54 

14.49 14.47 

14.43 

0.03 

0 

3 

SEDS J0333 16.95— 275338.8 

53.32061 -27.89410 

14.51 

14.44 

14.41 

14.38 

14.32 14.29 

14.23 

0.03 

0 

3 



14.86 

14.82 

14.69 

14.61 

14.57 14.56 

14.53 

0.03 

165 

1 

SEDS J033248. 1 1-273259.6 

53.20047 -27.54988 

14.65 

14.62 

14.55 

14.51 

14.50 14.49 

14.47 

0.03 

123 

1 



14.99 

14.98 

14.91 

14.87 

14.86 14.85 

14.83 

0.03 

205 

0 

SEDS J033219. 13— 273933.6 

53.07972 -27.65933 

14.68 

14.65 

14.60 

14.57 

14.56 14.56 

14.55 

0.03 

195 

1 



14.93 

14.92 

14.90 

14.89 

14.89 14.89 

14.88 

0.03 

411 

0 

SEDS J0333 15. 93— 280332.7 

53.31636-28.05909 

14.70 

14.67 

14.63 

14.59 

14.57 14.56 

14.54 

0.03 

0 

3 



14.77 

14.75 

14.74 

14.73 

14.72 14.72 

14.71 

0.03 

234 

1 

SEDS J0333 12. 34— 274232.8 

53.30142-27.70911 

14.73 

14.71 

14.69 

14.68 

14.67 14.67 

14.67 

0.03 

326 

1 



15.09 

15.09 

15.09 

15.09 

15.09 15.08 

15.08 

0.03 

405 

0 

SEDS J0333 18.60-2742 18.5 

53.32750 -27.70513 

14.73 

14.71 

14.68 

14.67 

14.67 14.66 

14.66 

0.03 

205 

1 



14.99 

14.99 

15.01 

15.01 

15.01 15.01 

15.01 

0.03 

335 

0 

SEDS J033 149. 13 — 2731 10.7 

52.95469 -27.51964 

14.75 

14.71 

14.61 

14.56 

14.54 14.52 

14.50 

0.03 

0 

3 



15.17 

15.11 

15.05 

15.01 

14.99 14.99 

14.97 

0.03 

0 

2 


Notes. As Table 7, but for the full-depth SEDS observations of the ECDFS. 

3 The PSF-fitted magnitude is given first, and the magnitudes given after are measured in apertures of 2"4, 3'.' 6, 478, 670, 772, and 1270 diameter. 
b Stated errors are la uncertainties for the 274 diameter aperture magnitudes. 

c Depth of coverage is expressed in terms of the number of unsaturated 100 s 3.6 fim IRAC frames that observed the source. 
d 3.6 fim photometric quality flag, as described in Section 5. If the flag is nonzero, the photometry may be corrupted. 
e Depth of coverage expressed in terms of the number of unsaturated 100 s 4.5 fim IRAC frames that observed the source. 
f 4.5 fim photometric quality flag, as described in Section 5. If the flag is nonzero, the photometry may be corrupted. 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.) 


Table 9 

SEDS Full-depth Source Catalog for the COSMOS Field 


Object 

R.A., Decl. 
(J2000) 


3.6 fim Magnitudes 3 
4.5 fim Magnitudes 


3.6 fim Unc. b 
4.5 fim Unc. 

3.6 /tm Cvg c 
4.5 Atm Cvg e 

3.6 fim Flag d 
4.5 fim Flag f 

SEDS J100032. 58+020825. 5 

150.13573 2.14043 

14.23 

14.06 

13.82 

13.64 

13.52 13.48 

13.38 

0.03 

69 

3 



14.09 

14.00 

13.76 

13.62 

13.56 13.54 

13.50 

0.03 

85 

3 

SEDS J100020. 87+01464 1.4 

150.08696 1.77816 

14.36 

14.29 

14.13 

14.03 

13.98 13.96 

13.92 

0.03 

35 

3 



14.61 

14.57 

14.42 

14.34 

14.30 14.29 

14.26 

0.03 

119 

1 

SEDS J100036.89+022357.4 

150.15371 2.39927 

14.39 

14.32 

14.17 

14.08 

14.03 14.01 

13.98 

0.03 

37 

3 



14.72 

14.67 

14.54 

14.46 

14.41 14.40 

14.36 

0.03 

80 

1 

SEDS J100017. 19+022554.8 

150.07163 2.43190 

14.40 

14.35 

14.23 

14.16 

14.12 14.11 

14.08 

0.03 

20 

3 



14.70 

14.67 

14.57 

14.51 

14.48 14.47 

14.45 

0.03 

59 

1 

SEDS J100032. 73+020505. 7 

150.13638 2.08493 

14.41 

14.34 

14.17 

14.07 

14.03 14.01 

13.97 

0.03 

27 

3 



14.72 

14.68 

14.52 

14.42 

14.38 14.36 

14.33 

0.03 

81 

1 

SEDS J 100027.70+022752.2 

150.11540 2.46451 

14.44 

14.39 

14.20 

14.12 

14.08 14.07 

14.03 

0.03 

49 

3 



14.76 

14.72 

14.57 

14.49 

14.45 14.44 

14.42 

0.03 

109 

1 

SEDS J100045. 09+021636. 8 

150.18787 2.27690 

14.50 

14.45 

14.31 

14.23 

14.19 14.18 

14.15 

0.03 

38 

1 



14.77 

14.74 

14.66 

14.60 

14.57 14.56 

14.54 

0.03 

72 

1 

SEDS J100041. 9 1+015912. 7 

150.17463 1.98685 

14.55 

14.51 

14.36 

14.29 

14.25 14.24 

14.21 

0.03 

116 

1 



14.90 

14.88 

14.78 

14.72 

14.69 14.68 

14.66 

0.03 

148 

0 

SEDS J100054.39+022723.9 

150.22663 2.45664 

14.61 

14.58 

14.47 

14.41 

14.39 14.38 

14.35 

0.03 

101 

1 



15.20 

15.15 

15.04 

14.95 

14.91 14.90 

14.87 

0.03 

23 

2 

SEDS J 100038.22+014730.0 

150.15927 1.79166 

14.67 

14.64 

14.54 

14.50 

14.47 14.46 

14.44 

0.03 

145 

1 



14.80 

14.80 

14.79 

14.77 

14.76 14.76 

14.76 

0.03 

263 

1 


Notes. As Table 7, but for the full-depth SEDS observations of the COSMOS field. 

a The PSF-fitted magnitude is given first, and the magnitudes given after are measured in apertures of 274, 376, 478, 670, 772, and 1270 diameter. 
b Stated errors are la uncertainties for the 274 diameter aperture magnitudes. 

c Depth of coverage is expressed in terms of the number of unsaturated 100 s 3.6 fim IRAC frames that observed the source. 
d 3.6 jLim photometric quality flag, as described in Section 5. If the flag is nonzero, the photometry may be corrupted. 
e Depth of coverage expressed in terms of the number of unsaturated 100 s 4.5 fim IRAC frames that observed the source. 
f 4.5 fim photometric quality flag, as described in Section 5. If the flag is nonzero, the photometry may be corrupted. 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.) 
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Table 10 

SEDS Full-depth Catalog for the HDFN 


Object 

R.A., Decl. 
(J2000) 


3.6 fim Magnitudes 3 
4.5 fim Magnitudes 


3.6 fim Unc b 
4.5 fim Unc. 

3.6 (im Cvg c 
4.5 (im Cvg e 

3.6 fim Flag d 
4.5 fim Flag f 

SEDS J123508.40+615934.9 

188.78499 61.99303 

14.06 

14.02 

13.91 

13.83 

13.76 13.72 

13.64 

0.03 

4 

3 



14.38 

14.28 

14.06 

13.95 

13.89 13.87 

13.82 

0.03 

1 

3 

SEDS J123737. 95+621630.4 

189.40813 62.27511 

14.13 

14.13 

14.08 

14.07 

14.07 14.06 

13.99 

0.03 

4 

3 



14.63 

14.55 

14.43 

14.37 

14.31 14.27 

14.20 

0.03 

4 

3 

SEDS J1235 13.78+620640.0 

188.80742 62.11110 

14.20 

14.18 

14.07 

13.99 

13.93 13.91 

13.86 

0.03 

23 

3 



14.66 

14.58 

14.42 

14.32 

14.26 14.24 

14.21 

0.03 

23 

3 

SEDS J123652. 97+620727.0 

189.22072 62.12416 

14.23 

14.21 

14.07 

13.97 

13.91 13.88 

13.80 

0.03 

0 

3 



14.52 

14.43 

14.25 

14.15 

14.10 14.08 

14.03 

0.03 

0 

3 

SEDS J123625. 04+6221 15.7 

189.10433 62.35435 

14.26 

14.25 

14.16 

14.10 

14.07 14.05 

14.03 

0.03 

84 

3 



14.82 

14.77 

14.65 

14.57 

14.53 14.52 

14.49 

0.03 

84 

1 

SEDS J123728. 27+622450.7 

189.36779 62.41407 

14.31 

14.29 

14.19 

14.13 

14.10 14.09 

14.06 

0.03 

232 

3 



14.83 

14.78 

14.64 

14.56 

14.52 14.51 

14.48 

0.03 

232 

1 

SEDS J123550.75+622608.9 

188.96144 62.43580 

14.39 

14.39 

14.33 

14.29 

14.27 14.26 

14.24 

0.03 

127 

3 



14.95 

14.89 

14.77 

14.68 

14.63 14.62 

14.59 

0.03 

127 

0 

SEDS J123845. 57+62045 1.2 

189.68988 62.08089 

14.45 

14.44 

14.36 

14.31 

14.29 14.28 

14.26 

0.03 

112 

3 



15.07 

15.02 

14.91 

14.83 

14.79 14.78 

14.76 

0.03 

56 

0 

SEDS J123633. 32+622400.9 

189.13884 62.40026 

14.52 

14.52 

14.43 

14.39 

14.37 14.36 

14.34 

0.03 

155 

1 



14.82 

14.79 

14.73 

14.69 

14.68 14.68 

14.67 

0.03 

318 

1 

SEDS J12382 1.37+621905.0 

189.58906 62.31806 

14.54 

14.34 

14.28 

14.21 

14.07 13.96 

13.56 

0.03 

43 

3 



14.73 

14.67 

14.60 

14.59 

14.58 14.52 

14.35 

0.03 

43 

3 


Notes. As Table 7, but for the full-depth SEDS observations of the HDFN. 

3 The PSF-fitted magnitude is given first, and the magnitudes given after are measured in apertures of 274, 376, 478, 670, 772, and 1270 diameter. 
b Stated errors are la uncertainties for the 274 diameter aperture magnitudes. 

c Depth of coverage is expressed in terms of the number of unsaturated 100 s 3.6 fim IRAC frames that observed the source. 
d 3.6 fim photometric quality flag, as described in Section 5. If the flag is nonzero, the photometry may be corrupted. 
e Depth of coverage expressed in terms of the number of unsaturated 100 s 4.5 fim IRAC frames that observed the source. 
f 4.5 fim photometric quality flag, as described in Section 5. If the flag is nonzero, the photometry may be corrupted. 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.) 


Table 11 

SEDS Full-Depth Source Catalog for the EGS 


Object 

R.A., Decl. 
(J2000) 


3.6 fim Magnitudes 3 
4.5 fim Magnitudes 


3.6 fim Unc b 
4.5 /im Unc. 

3.6 /tm Cvg c 
4.5 fi m Cvg e 

3.6 fim Flag d 
4.5 fim Flag f 

SEDS J1421 1 7.65+53095 1 .4 

215.32353 53.16427 

14.48 

14.47 

14.40 

14.35 

14.33 14.31 

14.29 

0.03 

42 

i 



15.04 

15.00 

14.89 

14.83 

14.79 14.78 

14.76 

0.03 

101 

i 

SEDS J141 83 1.98+523847.1 

214.63327 52.64643 

14.67 

14.66 

14.61 

14.57 

14.55 14.54 

14.52 

0.03 

127 

i 



14.93 

14.92 

14.89 

14.87 

14.86 14.85 

14.85 

0.03 

198 

i 

SEDS J141943. 89+524618.4 

214.93288 52.77177 

14.69 

14.69 

14.66 

14.64 

14.62 14.61 

14.60 

0.03 

94 

i 



14.99 

15.00 

15.00 

14.99 

14.98 14.97 

14.97 

0.03 

145 

i 

SEDS J142 102. 99+53 1520.5 

215.26245 53.25569 

14.70 

14.70 

14.65 

14.61 

14.59 14.58 

14.57 

0.03 

64 

i 



14.96 

14.96 

14.96 

14.95 

14.94 14.94 

14.94 

0.03 

103 

i 

SEDS J14 1727. 59+522043. 9 

214.36495 52.34553 

14.71 

14.70 

14.69 

14.68 

14.67 14.66 

14.65 

0.03 

0 

3 



14.88 

14.89 

14.90 

14.91 

14.91 14.92 

14.92 

0.03 

91 

1 

SEDS J141 846.52+524522.2 

214.69385 52.75618 

14.71 

14.71 

14.69 

14.68 

14.68 14.68 

14.67 

0.03 

201 

1 



15.05 

15.02 

15.04 

15.04 

15.04 15.04 

15.04 

0.03 

320 

1 

SEDS J14 1750. 60+523914. 4 

214.46085 52.65400 

14.72 

14.72 

14.70 

14.69 

14.68 14.68 

14.67 

0.03 

276 

1 



15.04 

15.04 

15.06 

15.07 

15.07 15.07 

15.07 

0.03 

337 

1 

SEDS J142 153. 23+53 1454. 6 

215.47180 53.24850 

14.79 

14.79 

14.80 

14.80 

14.80 14.80 

14.80 

0.03 

169 

1 



15.14 

15.15 

15.19 

15.20 

15.21 15.21 

15.23 

0.03 

163 

1 

SEDS J142107. 32+530030.0 

215.28052 53.00832 

14.81 

14.81 

14.75 

14.72 

14.70 14.69 

14.67 

0.03 

61 

1 



14.98 

14.98 

14.95 

14.94 

14.93 14.93 

14.92 

0.03 

115 

1 

SEDS J142101. 89+5303 16.3 

215.25786 53.05454 

14.81 

14.80 

14.79 

14.78 

14.78 14.78 

14.77 

0.03 

248 

1 



15.06 

15.07 

15.09 

15.10 

15.11 15.11 

15.12 

0.03 

413 

1 


Notes. As Table 7, but for the full-depth SEDS observations of the EGS. 

a The PSF-fitted magnitude is given first, and the magnitudes given after are measured in apertures of 274, 376, 478, 670, 772, and 1270 diameter. 
b Stated errors are la uncertainties for the 274 diameter aperture magnitudes. 

c Depth of coverage is expressed in terms of the number of unsaturated 100 s 3.6 fim IRAC frames that observed the source. 
d 3.6 jLim photometric quality flag, as described in Section 5. If the flag is nonzero, the photometry may be corrupted. 
e Depth of coverage expressed in terms of the number of unsaturated 100 s 4.5 fim IRAC frames that observed the source. 
f 4.5 fim photometric quality flag, as described in Section 5. If the flag is nonzero, the photometry may be corrupted. 

(This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.) 
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SEDS Magnitude (AB) 

Figure 29. A comparison of SEDS 2'.'4 diameter aperture photometry to that of four shallower coextensive surveys from the cryogenic mission. Symbols indicate the 
mean differences measured within 0.5 mag bins. All error bars are 1 cr . Panels (a) and (b): SEDS compared to SExtractor-based aperture photometry from SpUDS 
DR2. Vertical dotted lines indicate the SpUDS 5 a sensitivity limits. Panels (c) and (d): SEDS compared to SExtractor-based aperture photometry from S-COSMOS 
(Sanders et al. 2007). Vertical dotted lines indicate the 1 fj . Jy threshold of the S-COSMOS 3.6 // in selection, and the 1.7 //.ly S-COSMOS sensitivity at 4.5 fim. 
Panels (e) and (f): SEDS compared to the SExtractor-based aperture photometry from Barmby et al. (2008). The vertical dotted lines indicate the Barmby et al. (2008) 
50% completeness limits. Panels (g) and (h): SEDS compared to the SExtractor-based aperture photometry from SIMPLE version DR1. The vertical dotted lines 
indicate the 5cr point-source depths for the deepest SIMPLE integrations (2.5 hr). 


not coincide perfectly. Irregular regions on the periphery 
of the SEDS fields are covered deeply by only one of 
the two arrays. This version of the SEDS catalogs only 
includes sources imaged deeply in both of the SEDS 
passbands. 

5. The aperture photometry is likely to be corrupted for 
sources that fall on a strong negative residual from a 
relatively bright nearby source, and/or for sources that lie 
on or near pixels with (spurious) integration times of zero. 
Sources with photometry that may be corrupted owing to 
their proximity to at least one bright source have flags set 
to 1 in the catalogs. Sources with photometry that may be 
corrupted because of saturation have flags set to 2 in the 
catalogs. SEDS flags are determined separately for the 3.6 
and 4.5 /urn mosaics, and are cumulative. 

6. SEDS has been optimized specifically for faint-source 
photometry. Bright stars are problematic because they 
saturate in the 100 s SEDS exposures. Furthermore, SEDS 
photometry relies on a PSF-htting technique that cannot 
be expected to be accurate for extended sources (nearby 
galaxies). Extended sources will have more flux in the 
larger apertures than in the smaller ones. It is recommended 
that users rely on the 2'.'4 diameter aperture photometry for 
sources fainter than 21 mag because this aperture size is an 
appropriate match to the apparent solid angle of the faint 
sources that SEDS is designed to survey, and because this 
is also the aperture that best excludes contributions from 
nearby but unrelated sources in the field. 


6. DISCUSSION 
6.1. Confusion Noise 

The SEDS uncertainties are dominated by the 3% uncertainty 
in the IRAC absolute calibration for sources brighter than 
17 mag. They are likewise dominated by the IRAC sensitivity at 
the faint limit. At intermediate magnitudes, however, confusion 
noise is the dominant contributor to the uncertainties (Table 4). 
This is expected. The SEDS IRAC number counts meet the 
Rowan-Robinson (2001) criterion for confusion (1 source per 
40 beams) at roughly 21.8 mag in both bands. Confusion 
contributes to the errors even at brighter magnitudes, however. 
Between 18 and 25 mag, then, the uncertainties plotted in 
Figure 27 broadly reflect the confusion noise. Future surveys 
with highly redundant coverage may expect to encounter similar 
levels of confusion noise if their observations have the same 
angular resolution as SEDS. 

6.2. IRAC Color Distribution 

The IRAC colors of the sources give clues to their redshifts 
and luminosities. Sorba & Sawicki (2010) and Barro et al. 
(201 1) showed that the [3.6] — [4.5] color is a useful photometric 
redshift indicator, especially for separating galaxies at z < 1.3 
from those at z > 1.5. Figure 30 shows a color-magnitude 
diagram of the SEDS sources. As expected, sources with 
[3.6] < 19 have [3.6] — [4.5] ~ —0.4, consistent with late- 
type stars or low-redshift galaxies. Franceschini et al. (2006) 
found that galaxies in this magnitude range mostly have redshifts 
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Figure 30. Color-magnitude diagram for SEDS sources. The inset shows 
a diagram for sources with 14 < [3.6] < 18. For clarity, the main dia- 
gram only shows 1 source in 100 and the inset shows 1 in 10. There are 
302 objects (~0.1% of the catalog) with magnitudes and colors that put them 
outside the boundaries of both the main plot and the inset. 

between 0.5 and 1.2. There are many fainter objects with similar 
color, but starting around [3.6] 21, another population with 

[3.6] — [4.5] 0.1 appears. This is consistent with galaxies at 

z > 1.3 (Sorba & Sawicki 2010). Helgason et al. (2012) derived 
galaxy luminosity functions (LFs) at a range of wavelengths 
as a function of redshift. At z ~ 1.3, the turnover in the rest 
1.63 jxm LF, approximately corresponding to observed 3.6 /rm, 
is L* —24.0. This corresponds to apparent [3.6] ss 21.7, 
roughly a magnitude fainter than the brightest z > 1.3 galaxies. 
In other words, the brightest observed galaxies with colors 
consistent with z > 1.3 are roughly a magnitude brighter 
than L* at that redshift. The color trends are also shown by 
the histograms in Figure 31. At [3.6] > 23.5, precise color 
information is largely lost because of photometric uncertainties, 
but the histogram peak at [3.6] — [4.5] ~ —0.1 shows that the 
sources are a mix of low-redshift, low-luminosity galaxies and 
z > 1.3 galaxies with higher luminosities. 

6.3. Number Counts 

SEDS detects a total of roughly 300,000 sources in the 
1 .46 deg 2 covered by the survey. Figures 32 and 33 and Table 12 
present the resulting differential source counts along with Milky 
Way star counts estimated from the model of Arendt et al. 
(1998) for the SEDS lines-of-sight. At levels brighter than 
roughly 18 AB mag in both SEDS bands, the IRAC counts are 
consistent with the star count models. SEDS contains relatively 
few galaxies brighter than 18 mag. The vast majority of sources 
fainter than 18 mag, however, are galaxies: the contributions of 
Milky Way stars to the faint counts are negligible. 

As shown in Figures 32(e) and 33(e), the SEDS counts 
agree very well with those of Fazio et al. (2004a) down to the 
sensitivity limit of the Barmby et al. (2008) EGS survey, despite 
being based on a different detection/photometry technique. 
They are also consistent with the deeper counts in the HDFN 
(Figures 32(d) and 33(d)) compiled by Magdis et al. (2008). 



Figure 31. Histograms of SEDS source colors in four magnitude bins. The 
unlabeled histogram with the smallest number of sources is for [3.6] < 20. The 
inset shows the color distribution of red sources ([3.6] — [4.5] ^ 1). There are 
61 sources with [3.6] — [4.5] > 2 not shown in any histogram. All but one of 
these have [3.6] > 25.2. Magnitudes and therefore colors for sources that faint 
are highly uncertain. 

The field-to-field variation in number counts for objects with 
21 < [3.6] < 24 is only 3.3% of the area- weighted mean 
count (Figures 34(a) and (b), and Table 12). This has significant 
implications for the nature of the source population. Based on 
tests with the QUICKCV code (Newman & Davis 2002), such a 
low sample variance (or “cosmic variance” as it is commonly 
referred to) can only be achieved in one of two ways. One 
possibility is that the clustering of the sources is intrinsically 
weak (with a correlation length r (t ~ 3 h~ 1 Mpc comoving, 
matching low-mass star-forming galaxies at z = 0-1; Zehavi 
et al. 2011; Coil et al. 2008) and they have a moderately 
broad redshift distribution spanning A z ~ 1.5-2. 5. The other 
possibility is that the sources exhibit strong clustering, e.g., 
correlation length ro ~ 5-6 h~ 1 Mpc comoving, matching 
quiescent galaxies at z = 0-1 (Zehavi et al. 2011; Coil et al. 
2008) and Lyman-break-selected samples at z ~ 4 (Ouchi 
et al. 2004; Hildebrandt et al. 2009), but that collectively these 
strongly clustered sources are detected over a very wide redshift 
range, from z ~ 0-4. Subpopulations (e.g., color-selected to 
be at high redshift) show stronger cosmic variance than this, as 
would be expected if they have stronger clustering and/or span 
a smaller redshift range. 

If a single field of the same total area as SEDS were observed, 
errors in the mean number density or number counts of a 
population due to sample/cosmic variance would be a factor 
of ~ 1 .25 times higher than the errors that achieved for SEDS. 
SEDS yields lower variance because a single-field survey 
observes a correlated volume, which may be dominated by 


36 This number was obtained by comparing the predicted errors from QUICKCV 
for a single square field to the predicted error in the inverse-variance-weighted 
(i.e., optimal/minimum- variance) mean number density resulting from a set of 
five fields having the actual SEDS geometries. The actual improvement 
depends on the total redshift range spanned by a given SEDS sample, with 
single-field errors ranging from 1.21 times higher than SEDS for a sample 
exclusively at high z, 3.9 < z < 4.1; to 1.23 times higher if the objects span 
0 < z < 4; to 1.39 times higher for a lower redshift range, 0 < z < 2. 
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Mag 

(AB) 

UDS 


ECDFS 

COSMOS 

HDFN 


EGS 


Total 


Counts 

Unc. 

Counts 

Unc. 

Counts 

Unc. 

Counts 

Unc. 

Counts 

Unc. 

Counts 

Unc. 

3.6 fim 

14.75 

2.029 

0.243 

2.058 

0.224 

2.278 

0.236 

2.301 

0.200 

1.711 

0.333 

2.09 

0.32 

15.25 

2.077 

0.229 

2.012 

0.236 

2.136 

0.277 

1.806 

0.354 

2.012 

0.236 

2.02 

0.35 

15.75 

2.246 

0.189 

2.188 

0.192 

2.499 

0.183 

2.079 

0.258 

2.137 

0.204 

2.23 

0.24 

16.25 

2.317 

0.174 

2.410 

0.149 

2.420 

0.200 

2.182 

0.229 

2.380 

0.154 

2.35 

0.19 

16.75 

2.584 

0.128 

2.563 

0.125 

2.635 

0.156 

2.505 

0.158 

2.497 

0.135 

2.44 

0.14 

17.25 

2.547 

0.134 

2.665 

0.111 

2.800 

0.129 

2.471 

0.164 

2.596 

0.120 

2.61 

0.13 

17.75 

2.723 

0.109 

2.864 

0.088 

2.925 

0.112 

2.812 

0.111 

2.726 

0.104 

2.81 

0.10 

18.25 

3.061 

0.074 

3.128 

0.065 

3.146 

0.087 

3.052 

0.084 

3.034 

0.073 

3.08 

0.07 

18.75 

3.428 

0.048 

3.452 

0.045 

3.543 

0.055 

3.408 

0.056 

3.487 

0.043 

3.46 

0.04 

19.25 

3.775 

0.039 

3.726 

0.054 

3.839 

0.059 

3.791 

0.044 

3.802 

0.073 

3.78 

0.05 

19.75 

3.974 

0.031 

3.980 

0.035 

4.052 

0.041 

4.014 

0.032 

4.024 

0.044 

4.01 

0.03 

20.25 

4.130 

0.027 

4.142 

0.033 

4.223 

0.037 

4.149 

0.029 

4.175 

0.034 

4.16 

0.03 

20.75 

4.251 

0.025 

4.251 

0.032 

4.319 

0.037 

4.269 

0.026 

4.291 

0.032 

4.27 

0.03 

21.25 

4.357 

0.024 

4.360 

0.031 

4.411 

0.036 

4.392 

0.024 

4.409 

0.031 

4.38 

0.03 

21.75 

4.476 

0.023 

4.466 

0.031 

4.502 

0.035 

4.491 

0.024 

4.514 

0.029 

4.49 

0.02 

22.25 

4.584 

0.023 

4.579 

0.031 

4.599 

0.035 

4.582 

0.023 

4.606 

0.026 

4.59 

0.02 

22.75 

4.700 

0.023 

4.685 

0.021 

4.695 

0.022 

4.713 

0.023 

4.719 

0.026 

4.70 

0.02 

23.25 

4.813 

0.019 

4.787 

0.015 

4.800 

0.022 

4.816 

0.018 

4.818 

0.016 

4.81 

0.01 

23.75 

4.892 

0.017 

4.883 

0.015 

4.880 

0.023 

4.908 

0.019 

4.896 

0.016 

4.89 

0.01 

24.25 

4.949 

0.019 

4.913 

0.016 

4.916 

0.026 

4.994 

0.021 

4.917 

0.018 

4.94 

0.02 

24.75 

4.966 

0.048 

4.939 

0.047 

4.961 

0.040 

5.021 

0.036 

4.942 

0.052 

4.96 

0.04 

25.25 

5.022 

0.073 

4.951 

0.062 

4.966 

0.058 

4.979 

0.047 

5.029 

0.077 

4.99 

0.06 

25.75 

5.115 

0.169 

5.031 

0.108 

5.206 

0.131 

5.025 

0.078 

5.110 

0.152 

5.09 

0.14 

4.5 fim 

14.75 

1.799 

0.302 

1.602 

0.378 

2.136 

0.277 

1.903 

0.289 

1.058 

0.707 

1.76 

0.50 

15.25 

1.975 

0.258 

1.871 

0.277 

1.977 

0.333 

2.134 

0.243 

1.660 

0.354 

1.93 

0.46 

15.75 

2.160 

0.209 

2.079 

0.218 

2.136 

0.277 

1.944 

0.302 

2.099 

0.213 

2.09 

0.30 

16.25 

2.230 

0.192 

2.204 

0.189 

2.527 

0.177 

2.079 

0.258 

2.219 

0.186 

2.25 

0.23 

16.75 

2.330 

0.171 

2.390 

0.152 

2.485 

0.186 

2.283 

0.204 

2.400 

0.151 

2.38 

0.18 

17.25 

2.605 

0.125 

2.614 

0.118 

2.675 

0.149 

2.505 

0.158 

2.489 

0.136 

2.58 

0.14 

17.75 

2.674 

0.115 

2.721 

0.104 

2.897 

0.115 

2.556 

0.149 

2.716 

0.105 

2.71 

0.11 

18.25 

2.844 

0.095 

3.002 

0.075 

3.090 

0.092 

2.899 

0.101 

2.843 

0.091 

2.93 

0.08 

18.75 

3.265 

0.058 

3.238 

0.057 

3.330 

0.070 

3.240 

0.068 

3.269 

0.055 

3.27 

0.05 

19.25 

3.564 

0.075 

3.600 

0.050 

3.737 

0.063 

3.603 

0.081 

3.647 

0.050 

3.63 

0.06 

19.75 

4.069 

0.046 

3.894 

0.033 

3.993 

0.042 

3.946 

0.053 

3.979 

0.032 

3.98 

0.04 

20.25 

4.271 

0.043 

4.136 

0.028 

4.229 

0.037 

4.131 

0.049 

4.179 

0.028 

4.19 

0.03 

20.75 

4.382 

0.042 

4.260 

0.027 

4.323 

0.036 

4.290 

0.048 

4.295 

0.028 

4.31 

0.03 

21.25 

4.481 

0.044 

4.362 

0.026 

4.412 

0.036 

4.369 

0.047 

4.408 

0.027 

4.41 

0.03 

21.75 

4.579 

0.043 

4.445 

0.026 

4.488 

0.035 

4.462 

0.048 

4.493 

0.026 

4.50 

0.03 

22.25 

4.685 

0.043 

4.543 

0.025 

4.591 

0.035 

4.552 

0.047 

4.574 

0.026 

4.59 

0.03 

22.75 

4.778 

0.044 

4.655 

0.015 

4.673 

0.023 

4.642 

0.049 

4.673 

0.018 

4.69 

0.03 

23.25 

4.894 

0.025 

4.751 

0.014 

4.764 

0.022 

4.768 

0.019 

4.777 

0.014 

4.80 

0.02 

23.75 

4.975 

0.019 

4.841 

0.014 

4.856 

0.023 

4.842 

0.018 

4.849 

0.014 

4.88 

0.01 

24.25 

5.044 

0.020 

4.902 

0.015 

4.912 

0.026 

4.932 

0.020 

4.892 

0.015 

4.94 

0.02 

24.75 

5.150 

0.047 

4.974 

0.046 

5.004 

0.039 

5.050 

0.028 

4.990 

0.052 

5.03 

0.04 

25.25 

5.267 

0.073 

5.045 

0.062 

5.097 

0.056 

5.111 

0.030 

5.139 

0.076 

5.14 

0.06 

25.75 

5.158 

0.170 

5.179 

0.107 

5.112 

0.132 

5.237 

0.041 

5.128 

0.156 

5.17 

0.13 


Notes. Differential number counts measured for SEDS measured in bins of width 0.5 mag, expressed in terms of log (IV) mag” 1 deg” 2 . Counts given as “Total” are 
area- weighted means derived from all five SEDS fields using the areas given in Table 1. All uncertainties are la. The errors given for individual fields reflect only 
\fN counting errors, but the uncertainties attributed to “Total” counts also take field-field variations into account. 


common, large-scale high- or low-density fluctuations, whereas 
multiple widely separated fields sample statistically independent 
regions of the universe. This ratio of errors is independent 
of linear galaxy bias, and hence of the nature of the galaxies 
observed. 

While the ratio of errors between SEDS and a single field 
covering the total SEDS area is nearly independent of sample 


properties, the absolute reduction in errors will depend strongly 
on the source population’s bias and redshift range. As an 
example, we adopt the bias predictions given in Table 4 
of Moster et al. (2010) and their assumed power spectrum 
amplitude (erg = 0.77). In that case, a sample of galaxies 
with 10 8 ' 5 Mq < M < 10 9 M q which spans 0 < z < 2 would 
have 1.25% scatter in counts/number densities due to cosmic 
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14 16 18 20 22 24 26 

AB Magnitude [3.6 /im] 


Figure 32. Differential counts in the five SEDS fields at 3.6 fim. Open symbols show the raw counts, while solid symbols indicate the counts corrected for 
incompleteness on the basis of simulated detections as described in Section 4.2. Error bars represent the Poisson statistics only. The solid lines in panels (d) and (e) 
indicate the incompleteness-corrected counts measured in the HDFN by Magdis et al. (2008) and in the EGS by Fazio et al. (2004a), respectively. The dotted lines in 
all panels are the expected counts arising from Milky Way stars along the lines of sight from the DIRBE Faint Source Model at 3.5 /im (Arendt et al. 1998; Wainscoat 
et al. 1992; Cohen 1993, 1994, 1995). 


variance in SEDS, but would have a scatter of 1.60% in a single 
field of the same area. A sample with M > 10 11 M q and 
falling in the range 3.9 < z < 4.1 would have 21.0% scatter 
from sample variance in SEDS, versus 25.4% in a single-field 
survey. To achieve the same sample/cosmic variance as SEDS, 
a single-field survey of the same depth would have to cover 
~2.4 deg 2 , as opposed to the 1.46 deg 2 surveyed in the five 
SEDS fields. Of course, for some populations Poisson variance 
(shot noise) will dominate rather than sample variance; in that 
case only area surveyed matters. In general, the net errors in 
counts or number densities for a given sample from a single- 
field survey must be between 0% (if Poisson-dominated) and 
~25% (if sample-variance-dominated) worse than if the SEDS 
strategy were used, given that both terms contribute to the overall 
error. Given that observing time is (to first order) determined by 
the area surveyed, a multiple-field strategy will always yield 
errors that are no worse than, and can be significantly better 
than, a single-field observation at constraining the abundances 
of galaxy populations. 

The mean differential counts, calculated as area-weighted 
means, are shown in Figures 34(c) and (d). In the aggregate, 
the counts brighter than 18 mag in both bands very closely 
follow the Arendt et al. (1998) DIRBE 3.5 and 4.9 /urn star 
count models for the Milky Way. Elsewhere, the SEDS counts 
correspond closely to the “default model” counts of Helgason 
et al. (2012), who reconstruct galaxy counts in the IRAC bands 


from known galaxy populations using an extensive library of 
multi-wavelength LF measurements extending to z ~ 4-5. 

We fitted a pure power law to the SEDS counts below 
the break in both bands. Specifically, we performed a j 1 
minimization fit to the differential counts from 20-26 mag 
in both SEDS bands. The fitted slopes are 0.17 ±0.05 and 
0.180 ± 0.004 at 3.6 and 4.5 jim, respectively, where the uncer- 
tainties are la. This suggests that the dominant contributions 
to the faint IRAC counts come from galaxies at z — 3-5, with 
an LF having a faint-end slope of roughly a = — 1 .5 (Helgason 
et al. 2012). Although this slope is slightly lower than seen for 
UV LFs at similar redshifts ( a ~ —1.7; Bouwens et al. 2012; 
Bouwens et al. 2007), it is consistent with the trends in faint- 
end slope documented in Helgason et al. (2012) and anchored 
at low redshifts by the 3.6 and 4.5 /urn LFs of Dai et al. (2009) 
and Babbedge et al. (2006). 

6.4. The Integrated Galaxy Light from IRAC Sources 

The CIB contains the sum of all light from distant galax- 
ies integrated over all of Cosmic time, and therefore is of pri- 
mary importance for understanding galaxy formation and evo- 
lution. Historically, precise measurement of the CIB has been 
complicated by the difficulties involved in modeling and sub- 
tracting the very substantial foreground contaminants (see, e.g., 
the review by Kashlinsky 2005). However, measurement of the 
total light from resolved sources, the integrated galaxy light 
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AB Magnitude [4.5 jtim] 


Figure 33. As Figure 32, but for the differential 4.5 jitm SEDS counts, and showing the DIRBE Faint Source Model counts for Milky Way stars at 4.9 /xm. 


(IGL), provides a strict and robust lower limit on the CIB and is 
largely unaffected by the foregrounds that complicate direct CIB 
measurements. Low-background space-based surveys, such as 
SEDS, have obvious promise for this application. 

Fazio et al. (2004a) made the first measurements of the IGL 
with IRAC: 5.4 and 3.5 nW m 2 sr” 1 at 3.6 and 4.5 /x m respec- 
tively, down to 23.5 mag. Subsequent analysis by Kashlinsky 
(2005) revised these estimates respectively to 5.3 ±1.0 and 
4.0 ± 0.8 nW m 2 sr” 1 . In other words, the cryogenic IRAC 
observations resolved somewhat less than half of the estimated 
CIB at these wavelengths (13.3 ± 2.8 nW in 2 sr” 1 ; Levenson 
et al. 2007). 

The Fazio et al. (2004a) bright counts are drawn from wide- 
held observations covering nearly 10 deg 2 , and are therefore 
unlikely to suffer significantly from cosmic variance. However, 
as Levenson & Wright (2008) point out, the faint Fazio et al. 
(2004a) counts are based on only two relatively narrow lines 
of sight. Levenson & Wright (2008) therefore measured the 
IRAC counts toward three widely separated fields (DIRBE 
"Dark Spots”) each covering a substantial area (~1 deg 2 ) to 
estimate faint IRAC counts relatively free of cosmic variance. 
On the basis of an extrapolation of their counts to faint levels, 
Levenson & Wright (2008) arrived at a larger estimate of the 
contribution of IRAC sources to the CIB: 9.0 +l ( j y nW m” 2 sr -1 . 

SEDS offers a new opportunity — by virtue of its deep obser- 
vations spread over five widely-separated fields — to reduce the 
uncertainties attending the contribution of faint IRAC sources to 


the CIB. SEDS is of course not optimal for blight-source count 
estimates. We therefore adopt the Fazio/Kashlinsky measure- 
ments for the counts brighter than 23.5 mag. Despite the fact 
that SEDS reaches 2.5 mag deeper than Fazio et al. (2004a), the 
incremental contribution of faint SEDS sources (23.5-26.0 mag) 
to the previous CIB estimate is small, amounting to 
only ~10% of the earlier estimates: 0.33 ±0.03 and 
0.36 ±0.03 nW m 2 sr -1 at 3.6 and 4.5 /xm, respectively. The 
total contributions to the CIB of detected IRAC sources are 
therefore 5.6 ± 1.0 and 4.4 ± 0.8 nW m” 2 sr -1 at 3.6 and 
4.5 /xm, respectively, down to 26 mag. Our results agree with 
Helgason et al. (2012), whose “default” construction estimated 
total contributions from IRAC sources to the CIB of 4 .Wq 7 7 and 
3-3 ^q 2 5 nW m” 2 sr" 1 at 3.6 and 4.5 /xm, respectively, where the 
ranges bracket the two extremes of the faint-end LF extrapola- 
tion consistent with the Helgason et al. (2012) LF database. 

The SEDS CIB measurement is significantly smaller than 
that estimated in the three DIRBE “Dark Spots” at 3.6 /xm by 
Levenson & Wright (2008). This may be because the “Dark 
Spots” estimate is based on an extrapolation to faint levels 
of relatively bright counts having a significantly steeper slope 
(~0.5) than found by SEDS. Figure 12 of Levenson & Wright 
(2008) nonetheless suggests that the measured SEDS slope is 
consistent with their analysis at roughly the 90% confidence 
level. 

The faint SEDS IRAC counts impose strong upper limits on 
the possible contribution of IRAC sources to the CIB. Even if 
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Figure 34. Differential source counts for SEDS. Panel (a): symbols indicate the 3.6 /zm counts in the five SEDS fields. The black lines indicate the counts after 
correcting for incompleteness. The red line shows the incompleteness-corrected counts from Fazio et al. (2004a). Panel (b): as panel (a), but for the 4.5 /zm counts. 
Panel (c): the mean differential 3.6 /zm source counts for all of SEDS, weighted by field area. The upper and lower dashed blue lines indicate respectively the 
high-faint-end (HFE) and low-faint-end (LFE) count models from Helgason et al. (2012). The dashed blue line between these two models indicates the default 
Helgason et al. (2012) model, which averages the HFE and LFE models. The black dotted line shows the area-weighted mean predicted Milky Way star counts based 
on the DIRBE 3.5 /zm Faint Source Model Counts (Arendt et al. 1998). Panel (d): as panel (c), but for the 4.5 /zm counts and showing the DIRBE Faint Source Model 
counts for Milky Way stars at 4.9 /zm. 

(A color version of this figure is available in the online journal.) 


the slopes of the IRAC counts increased immediately below 
26 AB mag to the Helgason high faint end case (0.22 in both 
SEDS bands), the total contributions to the CIB from these 
undetected sources would be only 0.16 and 0.14 nW m -2 sr~ 1 
at 3.6 and 4.5 fim, respectively. Thus, discrete IRAC sources 
appear unable to furnish more than about half of the estimated 
CIB at 3.6 and 4.5 /xm. 

6.5. Other Applications of SEDS Data 

SEDS is designed to maximize its impact in fields related 
to cosmology and high-redshift science. It has already been 
of some use in this regard. Kashlinsky et al. (2012) used 
SEDS imaging from two fields (the UDS and EGS) to make 
new measurements of fluctuations in the CIB; the SEDS data 
permitted extraction of much fainter sources and analysis of 
the CIB fluctuations over significantly larger scales than was 
possible previously at these wavelengths. 

SEDS will be of particular interest as a means to identify 
distant galaxies. High-redshift galaxies have long been identified 
on the basis of their red colors; historically this has been 
accomplished in the optical-to-NIR regime (e.g., I — K or R — K) 
to bracket galaxies’ 4000A break and measure their redshifts 
out to z ~ 3 (Roche et al. 2003; Franx et al. 2003; Brown et al. 
2005). With Spitzer / IRAC photometry this technique has been 
extended out beyond z = 3 (e.g., Wang et al. 2009; Huang 
et al. 2011; Wang et al. 2012). This has begun happening with 


SEDS photometry as well. Aravena et al. 2010 incorporated 
SEDS photometry into the spectral energy distribution of a rare 
submm-selected galaxy in COSMOS to place it at z ~ 4.8. 
Caputi et al. (2012) identified and analyzed a sample of very rare, 
very red objects ( H — [4.5] > 4) galaxies in the UDS and found 
that this color criterion was efficient for identifying sources at 
3 < z < 5. At an even more extreme redshift, Capak et al. 
(2011) used SEDS-COSMOS observations (along with a suite 
of other data) to identify a massive protocluster at z ~ 5.3. No 
doubt the SEDS data will be put to other uses not yet anticipated 
by the SEDS team, now that they are public. 

7. SUMMARY 

The SEDS team has carried out a deep 3.6 and 4.5 pm survey 
of five widely separated extragalactic fields, covering a total of 
1.46 deg 2 to a depth of 26 AB mag (3er) in both bands. The 
images and catalogs being made publicly available with this 
contribution therefore constitute a significant advance into a new 
regime of depth-area survey parameter space, never previously 
explored at these wavelengths. 

SEDS has detected and photometered more than 300,000 dis- 
tinct sources, of which the vast majority are galaxies. The galaxy 
color-magnitude distribution is consistent with mostly z < 1.3 
galaxies at [3.6] < 21 and a mix of low- and high-redshift 
galaxies at fainter magnitudes. Detailed number count models 
will have to reproduce both the number counts themselves and 
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the color distributions as a function of magnitude. The resulting 
counts are consistent with but extend previous measurements 
to significantly fainter levels, and are consistent with infrared 
galaxy LFs having faint-end slopes described by power laws 
with indices of a = —1.5. The contribution of resolved IRAC 
sources to the CIB is approximately half of the diffuse back- 
ground measured by DIRBE. 
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